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C<"> ■ Abstract 
> 

£ — , Asymptotic spectral distribution (ASD) of the crosscorrelation matrix is investigated for a random 

O ■ 

spreading short/long-code asynchronous direct sequence-code division multiple access (DS-CDMA) sys- 
tem. The discrete-time decision statistics are obtained as the output samples of a bank of symbol matched 
filters of all users. The crosscorrelation matrix is studied when the number of symbols transmitted by 
q , each user tends to infinity. Two levels of asynchronism are considered. One is symbol-asynchronous but 

chip-synchronous, and the other is chip-asynchronous. The existence of a nonrandom ASD is proved by 
moment convergence theorem, where the focus is on the derivation of asymptotic eigenvalue moments 
(AEM) of the crosscorrelation matrix. A combinatorics approach based on noncrossing partition of set 
partition theory is adopted for AEM computation. The spectral efficiency and the minimum mean-square- 
error (MMSE) achievable by a linear receiver of asynchronous CDMA are plotted by AEM using a 
numerical method. 



'This paper was presented in part at IEEE ISIT'05, Adelaide, Australia, Sept. 4-9, 2005 and IEEE ISIT07, Nice, France, 
June 24-29, 2007. 
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I. Introduction 

Direct sequence-code division multiple access (DS-CDMA) is one of the most flexible and commonly 
proposed multiple access techniques for wireless communication systems. To gain deeper insights into 
the performance of receivers in a CDMA system, much work has been devoted to the analysis of random 
spreading CDMA in the large-system regime, i.e., both the number of users K and the number of chips 
N per symbol approach infinity with their ratio K/N kept as a finite positive constant (3 [l]-[3]. Such 
asymptotic analysis of random spreading CDMA enables random matrix theory to enter communication 
and information theory. In the last few years, a considerable amount of CDMA research has made 
substantial use of results in random matrix theory (see [4] and references therein). 

In this paper, we focus on the derivation of the asymptotic spectral distribution (ASD) of the cross- 
correlation matrix in asynchronous CDMA systems. Consider the linear vector memoryless channels of 
the form y = Hx+w, where x, y and w are the input vector, output vector and additive white Gaussian 
noise (AWGN) vector, respectively, and H denotes the random channel matrix independent of w. This 
linear model encompasses a variety of applications in communications such as multiuser channels, multi- 
antenna channels, multipath channels, and, in particular, asynchronous CDMA channels of our interest 
in this research, etc., with x, y and H taking different meanings in each case. Concerned with the linear 
model, it is of particular interest to investigate the limiting distribution of eigenvalues of the random 
matrix HH^ or H^H, called the ASD of the random matrix, when the column and row sizes of H tend 
to infinity but the ratio of sizes is fixed as a finite constant. Since ASD is deterministic and irrelevant 
to realizations of random parameters, it is convenient to use the asymptotic limit as an approximation 
for finite-sized system design and analysis. Moreover, it is quite often that ASD provides us with much 
more insights than an empirical spectral distribution (ESD) does. Even though ASD is obtained with the 
large-system assumption, in practice, the system enjoys large-system properties for a moderate size of 
the channel matrix. 

Some applications of ASD in communication and information theory are exemplified below. Take the 
linear model y = Hx + w for illustration. A number of the system performance measures, e.g. channel 
capacity and the minimum mean-square-error (MMSE) achievable by a linear receiver, is determined 
by the ESD of the matrix HHK The asymptotic capacity and MMSE obtained by using ASD as an 
approximation of ESD can often result in closed-form expressions [2], [5]. It is also shown in [6]-[9] 
that empirical eigenvalue moments (or, more conveniently, AEM) can be used to find the optimal weights 
of the reduced-rank MMSE receivers and its output signal-to-interference -ratio (SIR) in a large system. 
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Moreover, a functional related to AEM is denned as the free expectation of random matrices in free 
probability theory [10], which has been recently applied to the asymptotic random matrix analysis. 

For synchronous DS-CDMA systems, it is well known that the ASD of the crosscorrelation matrix 
follows Marcenko-Pastur law [5]. Also, explicit expressions for the AEM under the environments of 
unfaded, frequency-flat fading with single and multiple antennas, and frequency-selective fading are 
derived in [11], [12]. Actually, most of the research results on random spreading CDMA making use of 
random matrix theory are confined to synchronous systems. Just a few of them investigate asynchronous 
systems [13]— [19]. The goal of this work is to find out the ASD of crosscorrelation matrix in asynchronous 
CDMA systems given a set of users' relative delays and an arbitrary chip waveform. As the uplink of a 
CDMA system is asynchronous, this work is motivated by the needs to study the problem of asynchronous 
transmission that is important but much less explored in the area of random matrix theory. 

Two levels of asynchronism are considered in this paper, i.e., symbol-asynchronous but chip-synchronous, 
and chip-asynchronous. In the sequel, chip-synchronous is used for short to denote the former, and symbol- 
synchronous represents an ideal synchronous system. To be more specific, the relative delays among users 
are integer multiples of the chip duration in chip-synchronous CDMA, while they are any real numbers 
in chip-asynchronous CDMA. 

Some previous results on asynchronous CDMA are reviewed. In [13], it is shown that the output SIR of 
the linear MMSE receiver in chip-synchronous CDMA converges to a deterministic limit characterized 
by the solution of a fixed-point equation that depends on the received power and the relative delay 
distributions of the users. When the width of the observation window during detection tends to infinity, the 
limiting output SIR converges to that of a symbol-synchronous system having all identical parameters. The 
system model of [13] splits each interferer into two virtual users, which leads to a crosscorrelation matrix 
with neither independent nor identically distributed entries. Results of [13] are obtained by employing 
Stieltjes transform for random matrices of that type. In [16], some equivalence results about the MMSE 
receiver output are provided for CDMA systems with various synchronism levels. In specific, when the 
ideal Nyquist sine pulse is adopted as the chip waveform of a chip-asynchronous system, the asymptotic 
SIR at the MMSE detector output is the same as that in an equivalent chip-synchronous system for any 
observation window width; moreover, as the observation window width increases, the output SIR in chip- 
asynchronous CDMA converges further to that in an equivalent symbol-synchronous system. In [18], the 
analysis of linear multiuser detectors is provided for a symbol quasi-synchronous but chip-asynchronous 
system, called symbol-quasi-synchronous for short. It is demonstrated that, when the bandwidth of the 
chip waveform is smaller than 1/(2T C ), where T c is the chip duration, the performance of symbol-quasi- 
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synchronous and symbol-synchronous systems coincides independently of the chip waveform and the 
distribution of relative delays among users, where the performance is characterized by the output SIR of 
a reduced-rank MMSE detector when a square-root raised cosine pulse is adopted as the chip waveform. 
If the bandwidth is larger than the threshold, the former system outperforms the latter. Actually, when the 
chip waveform bandwidth is narrower than the threshold, the inter-chip interference (ICI) free property 
is lost [20], which leads to a severe degradation in performance. 

In this paper, the system model is constructed without the user splitting executed in [13]. In stead, 
sufficient statistics obtained in the same way as [21]-[23] are employed. We assume the width of the 
observation window for symbol detection tends to infinity. The formulas for AEM of the crosscorrelation 
matrix are derived using a combinatorics approach. In specific, we use noncrossing partition in set partition 
theory as the solving tool to exploit all nonvanishing terms in the expressions of AEM. The combinatorics 
approach has been adopted in [1 1], [24]-[28] to compute AEM of random matrices in symbol-synchronous 
systems. All of them, either explicitly or implicitly, make use of graphs to signify noncrossing partitions. In 
this work, a graphical representation of K-graph, which is able to simultaneously represent a noncrossing 
partition and its Kreweras complementation map [29], is adopted. This property of a K-graph facilitates 
the employment of noncrossing partition and free probability theory in solving problems of interest. 

In some applications of probability theory, it is frequent that the (infinite) moment sequence of an 
unknown distribution is available, and these moments determine a unique distribution. Suppose that the 
goal is to calculate the expected value of a certain function g of the random variable X whose distribution 
is unknown. One of the most widely used techniques is based on the Gauss quadrature rule method [30], 
where the expected value of g(X) is expressed as a linear combination of samples of g(x), and moments 
of X are used to determine the coefficients in the combination and the points to be sampled. In this 
paper, the Gauss quadrature rule method is employed to compute the spectral efficiency and MMSE of 
asynchronous CDMA using the derived AEM. 

This paper is organized as follows. In Section II, the crosscorrelation matrices are given for chip- 
synchronous and chip-asynchronous CDMA systems. Some definitions regarding the limit of of a random 
matrix are also introduced. In Section III, we derive AEM and ASD of corsscorrelation matrices in both 
chip-synchronous and chip-asynchronous systems. In Section IV, free probability theory is employed to 
obtain the spectra of sum and product of crosscorrelation matrix and a random diagonal matrix. In Section 
V, mathematical results demonstrated in this paper are connected to some known results. Discussions of 
the spectral efficiency and MMSE in an asynchronous CDMA system are provided in Section VI. Finally, 
this paper is concluded in Section VII. 
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II. Crosscorrelation Matrix of Asynchronous CDMA 

Consider asynchronous direct sequence-code division multiple access (DS-CDMA) systems where each 
user's spreading sequence is chosen randomly and independently. There are K users in the system, and 
the number of chips in a symbol is equal to N. We focus on the uplink of the system and assume 
the receiver knows the spreading sequences and relative delays of all users. Systems with two levels of 
asynchronism are considered, i.e., symbol-asynchronous but chip-synchronous, and chip-asynchronous. 
In the sequel, chip-synchronous is used for short to denote the former, and symbol-synchronous refers 
to an ideal synchronous system. To differentiate notations of chip-synchronous and chip-asynchronous 
systems, subscripts in text form of "cs" and "ca" are used for notations in the former and the latter 
systems, respectively. 

A. Chip-Synchronous CDMA 

Denote the relative delay of user k as r^. For convenience, users are labelled chronologically by their 
arrival time, and {jk}^ =1 satisfy 

= Ti < t 2 < ■ ■ ■ < t k < NT C , (1) 

where T c is the chip duration, and all r^'s are integer multiples of T c . Suppose that each user sends a 
sequence of symbols with indices from — M to M. In the complex baseband notation, the contribution 
of user / to the received signal in a frequency-flat fading channel is 

M (n+l)N-l 
Xl {t)= MnMn) J2 ^^(t-qTc-n), 

n=-M q=nN 

where bi(n) is the n-th symbol of user I, Ai{n) is the complex amplitude at the time bi(n) is received, 
{of 1 : nN < q < (n + l)N — 1} is the spreading sequence assigned to the n-th symbol of user I, and 
ip(t) is the normalized chip waveform having the zero ICI condition of 



j — < 



f 1, r = 0, 

(2) 

0, r£Z\{0}. 



4>(t)iP(t - rT c )dt = 

J — oo 

It is assumed that {bi(n)}^L_ M is a collection of independent equiprobable ±1 random variables. The 
symbol streams of different users are independent. The distribution of a scaled chip y/Nc^ has zero 
mean, unit variance and finite higher order moments. We do not assume a particular distribution of 
cf \ Two distinct spreading mechanisms are considered, i.e., short-code and long-code. In a short-code 
system, the spreading sequences are randomly chosen for the first symbols, i.e., bk(—M) for user k, and 
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remain the same for every symbol. In a long-code system, the spreading sequences are randomly and 
independently picked from symbol to symbol. The sequences of received amplitudes {A k {m)}^ = _ M 
and {Ai(n)}^L_ M are independent if k ^ I. 

The complex baseband received signal is given by 

K 

r(t) = +«>(*)> (3) 

l=i 

where w(t) is the baseband complex Gaussian ambient noise with independent real and imaginary 
components. The correlation function of w(t) is E{w(t)w(t + r)*} = N 5d(t) with <5d(t) being the 
Dirac delta function. The symbol matched filter output of user fc's symbol m, denoted as yk(m), is 
obtained by correlating r(t) with the signature waveform of user fc's symbol m 

,oc /(m+l)JV-l \ 

Vk(m) = / r(t) Yl c { ^{t-pT c - Tk ) \ dt (4) 

J -°° \ P =mN J 

K M 

= E E Ai(n)bi(n)p cs (m,n;k,l) +v k (m), (5) 

1=1 n =-M 

where Vk(m) results from the ambient noise w(t), and p cs (m, n; k, I) is the crosscorrelation of spreading 
sequences at user fc's m-th symbol and user Vs n-th symbol, given as 

(n+l)JV-l (m+l)JV-l 

E E c f 9) 4 P) / - ^ - - P T c - r k )dt. (6) 

g =nAf p=mN ^ 00 

Due to the zero ICI condition of (2), the integral in (6) is nonzero (equal to one) if and only if pT c + Tk = 
qT c + t\. Thus, we obtain 

(m+l)JV-l (tH-I)JV-I 

p cs {m,n;k,l) = ^ ^ c^c^ 5( P T C + r k , qT c + n), (7) 

p=mN q=nN 

with 5(i,j) the Kronecker delta function. Since < Tfe,r; < (iV — 1)T C , for a specific symbol index m, 
the <5 function in (7) is equal to zero if n £ {m — 1, m, m + 1}. Thus, we rewrite (5) as 

K min{m+l,M} 

= E E Ai(n)k(n)p cs (m,n;k,l) + v k (m), —M < m < M. (8) 

2=1 n=max{m- 1,— M} 

Define the symbol matched filter output vector at the m-th symbol as 

y(m) = [yi(m),y 2 {m),--- ,y K (m)] T , 

and let the transmitted symbol vector b(m) and the noise vector v(m) have the same structures as 
y(m). Moreover, we define a block matrix R Cfi whose (k,l)-th element of the (m,n)-th block, with 
— M < m,n < M, 1 < k, I < K, is equal to p Cfi (m,n; k,l) in (7). The square bracket [•] is used to 
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indicate a specific element of the matrix R cs . Specifically, [R cs ]mn,ki represents the (k,l)-th entry of the 
(m, n)-th block of the block matrix ,R CS . When we just want to point out a specific block, only the first 
set of indices is used, i.e., [R cs ]mn- 

Using the notations defined above, we can show from (8) that 

min{m+l,M} 

y(m) = ^2 [Rcs]mnA(n)b{n) + v(m), 

ra=max{m— 1,— M} 

where A(n) = diagjA^n), A 2 (n), • • • , Ax{n)}. Stacking up y(m)'s to yield the symbol matched filter 
output of the whole transmission period as 

y=\y T (-M),y T (-M + l),--- ,y T (M)] T , 

we obtain the discrete-time signal model 

y = R cs Ab + v, (9) 

where b and v have the same structures as y, A = diag{A(— M), A{— M + 1) • • • , A(M)}, and the block 
matrix R Cfi has a tri-diagonal structure of 



[■Rcs]-l-2 [-Rcsl-1-1 [-RcsJ-10 

csJOl 

[-Rcs] 10 [-Rcs] 11 [-Rcs] 12 



(10) 



Since < t\ for k < I, [R Cii } m m -i and [^ C s]m m+i are strict (zero diagonal) upper- and lower-triangular 
matrices, respectively. From the signal model given in (9), .R cs can be viewed as the crosscorrelation matrix 
of chip-synchronous CDMA. It can be shown that the correlation matrix of the noise vector v in (9) is 
Elm/ 1 } = iVo-Rcs- Let ^ cs = H cs Hl s be a decomposition of ,R CS . We can perform the whitening process 
by left-multiplying y in (9) with if" 1 , resulting in 



z = Hjy = HlAb + w, 
where the noise vector w = H~^v is white and has the correlation matrix iVo-f- 



(11) 



B. Chip-Asynchronous CDMA 

In chip-asynchronous CDMA, the assumption that r^'s are integer multiples of the chip duration T c 
no longer exists. Although the relative delays of users are not integer multiples of T c , it is assumed that 
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the zero ICI condition still holds. Thus, the property of zero ICI exists for chips of each particular user. 
Similarly to (5), the symbol matched filter output yk(m) can be expressed as 

K M 

yk(m) = ^2 X A l(n)bi(n)p ca (m, n; k, I) + v k (m), (12) 

1=1 n=~M 

where p ca (m, n; k, I) is different from p cs (m, n; k, I) in (5) since the zero ICI condition does not hold when 
the time difference of chip waveforms is not integer multiples of T c . At this moment, the crosscorrelation 

/9 ca (m, n; k, I) is given by 

(m+l)AT-l (n+l)W-l 

p c ,(m,n;k,l) = ^ X c[ p) R^p - q)T c + r k - n) , (13) 

p=mN q=nN 

where 

/oo 
ip(t)ij(t - x)dt (14) 
-oo 

is the autocorrelation function of the chip waveform ip(t). We define the block matrix R ca whose (k, Z)-th 
component of the (m, n)-th block, with — M < m,n < M and 1 < k, I < K, is equal to p ca (m, n; k, I) 
given in (13). It can be shown that 

M 

y(m)= ^ [R Cii \ mn A(n)b(n) +v(m), (15) 

n=-M 

and we obtain the discrete-time signal model 

y = R CSL Ab + v, (16) 

where R ca is thus seen as the crosscorrelation matrix of a chip-asynchronous CDMA system. Note that, 
unlike the tri-diagonal structure of ,R CS shown in (10), the matrix R ca generally does not possess such 
structure except when the autocorrelation function R^(x) has a finite span. We can perform whitening 
on (16) to yield a linear model 

z = HlAb + w, (17) 

where il ca = H ca H\ a and w = H^v is a white noise vector. 

The discrete statistics y in (9) and (16), and hence z in (11) and (17), for chip-synchronous and 
chip-asynchronous systems, respectively, are sufficient and are obtained in the same way as [21]-[23]. 
These sufficient statistics are the output samples of a bank of filters matched to the symbol spreading 
waveforms of all users. An alternative approach to generating statistics, adopted by [13], [15], [16], 
[31], is to pass the received signal to a chip matched filter and sample the output. Statistics yielded in 
this way are sufficient only under symbol- and chip-synchronous assumptions, and are not sufficient in 
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the chip-asynchronous case. For a chip-asynchronous system, it is reported in [32] that, when the chip 
waveform is time limited to the chip interval, statistics obtained by sampling the chip matched filtering 
output at the chip rate leads to significant degradation in performance. On the other hand, if the output 
is sampled at up to the Nyquist rate, the loss in performance is negligible. 

C. Asymptotic Spectral Distribution of Crosscorrelation Matrix 

The analysis of asynchronous CDMA systems will be conducted in a large system regime. That is, we 
assume both the number of users K and the spreading gain N approach infinity with their ratio K/N 
converging to a non-negative constant j3. To proceed the analysis, some definitions regarding the limit 
of a random matrix [33] are introduced. Let denote a p x p Hermitian random matrix whose each 
element is a random variable. Suppose that has eigenvalues v\ < v-i < ■ ■ ■ < v p . Since S [p > is 
Hermitian, all z/j's are real. The ESD of is defined as 

F(P\x)=p- 1 #{i:v i <x}, (18) 

where #{•••} denotes the number of elements in the indicated set. A simple fact is the n-th moment of 
(x) can be represented as 

/oo 
x n AF ( P\x) = p- 1 tr((5 (p) ) n ), (19) 
-oo 

where tr is the trace operator, and the second equality holds because tr(( 1 S ,( - p - ) ) n ) = Yli=i • ^ F^ p \x) 
converges to a nonrandom function F(x) as p — > oo, then we say that the sequence {S^ : p = 1, 2, ■ • • } 
has an ASD F(x). To show that F^ p \x) tends to a limit, the moment convergence theorem [34] can be 
employed. To be specific, the theorem is stated here in a form convenient for this paper. 

Theorem 1: [Moment Convergence Theorem] Let {F^ p \x) : p = 1, 2, • ■ • } be a sequence of distribu- 
tion for which the moments 

/oo 
x n dF^(x) 
-oo 

exist for all n = 0, 1, 2, • • • . Furthermore, let F(x) be a distribution function for which the moments 

/oo 
x n dF{x) (20) 
-oo 

exist for all n = 0, 1, 2, • • • . If 

lim = m n (21) 

p^oo 
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for all n = 0, 1, 2, • • • in some sense, and if F(x) is uniquely determined by the sequence of moments 
mo, mi, 7U2, • • • , then 

lim F^\x) = F(x), 

and the convergence holds in the same sense as that of (21). ■ 
The moment convergence theorem has a long history. The details can be found in [35]. In applying this 
theorem to show the existence of the ASD F(x), we should determine the asymptotic moment sequence 
{m n } in (20) and prove that a unique distribution is determined by {m n }. In [36], Carleman gave a 
sufficient condition XmLi m 2n^ 2 "^ = °°> cai l e d Carleman's criterion, for the uniqueness of a distribution 
given a moment sequence {m n }. 

Concerned with the linear models of (11) and (17), it is of interest to consider the ESD of the random 
matrix H\A(H\A)\ x G {cs,ca} [4, Chapter 1]. Represent matrices ff x , A, and R x by H[ K \ A^ k \ and 
Ri K \ respectively, when the user size of the system is K. In order to use Theorem 1 to find the ASD 
of the random matrix it is required to find the limits of the empirical moments 

(2M + l)- 1 r 1 tr((ffW t #U( Jf ) t lf( Jf )) B ), 
= {2M + l)- 1 K- 1 tr{(A^R[ K '>A^ K '>) n ), n>0. (22) 

In an unfaded channel, i.e., \Ak(m)\ 2, s are equal for all k and m, the matrix A^ is a scaled identity 
matrix. Thus, the quantity 

(2M + l)- l K-hv{(R { p) n ) (23) 

is of interest in (22). 

III. ASD of Crosscorrelation Matrix in Asynchronous CDMA 

The goal of this section is to show that the ESD's of A^RiPA^ and A^R^A^ converge 
to nonrandom limits when K, N — > oo and K/N — > [3. We consider chip-synchronous and chip- 
asynchronous CDMA in Sections III-A and III-B, respectively. 

A. Chip-Synchronous CDMA 

We will use Theorem 1 to prove the result stated in the previous paragraph. We first consider the 
case of unfaded channel, and then we extend to a frequency-flat fading channel. The proof starts from 
showing the existence of 

= lim (2M + ly'K-'E {tr(GR(f >)»)} , n > 0, (24) 

K/N-+/3 
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where the functional fi(-) is a limiting normalized expected trace of the matrix in the argument, and the 
limit M — > oo is placed because we investigate the system behavior when the width of the observation 
window for symbol detection tends to infinity. We take the relative delays r^'s as deterministic quantities, 
and the expectation of (24) is with respect to (w.r.t.) the random spreading sequences. We have the 
following lemma. 

Lemma 1: In both short-code and long-code chip-synchronous CDMA systems, for any relative delays 
{Tk}k = i and any chip waveform ijj(t) satisfying (2), /i(^™s) exists and is given by 

(25) 

Proof: See Appendix II. ■ 
In Appendix II, we prove this lemma with the aid of techniques from noncrossing partition. Results 
of noncrossing partitions necessary for the proof are summarized in Appendix I. The same tool has 
been employed in [11,24] for a symbol-synchronous system. Note that, in the proof of Lemma 1, the 
spreading sequences are only assumed to be independent across users. For a particular user, we do not 
assume that the sequence is independent across symbols. Thus, the proof is applicable to both short-code 
and long-code systems. Moreover, the relative delays {jk}^ =1 are treated as deterministic constants, and 
we do not adopt a particular chip waveform except for the zero ICI condition. Thus, /i(i?c S ) does not 
depend on the asynchronous delays and the chosen chip waveform. 

Lemma 2: The n-th moment of the ESD of converges a.s. to ^(H™) when M, K, N — ► oo 

and K/N — > [3. Moreover, the moment sequence {/j,(R™ s ) : n > 1} satisfies the Carleman's criterion 

E~iMier l/(2n) =°o. 

Proof: See Appendix III. ■ 
Since the n-th moment of the ESD of Rffl converges to /j,(R™ s ), we refer to ^(R^) as the n-th 
AEM of R cs . It is seen that fi(R^ s ) given in (25) is equal to the n-th moment of the Marcenko-Pastur 
distribution [37] with ratio index (3, having density 

/,<*>= (i-±r*M + ^-' ,)+( >-* )+ 



(5) v ' 2ir(3x 
where (z) + = max{0,z}, a = (1 — \//3) 2 and b = (1 + \f^) 2 ■ As Lemma 2 shows the n-th empirical 
moment of R^ converges to n(Rj^) for n > and the moment sequence {/z(i?" s )} satisfies Carleman's 
condition, the following theorem holds straightforwardly due to Theorem 1. 

Theorem 2: In both short-code and long-code chip-synchronous CDMA systems, for any relative delays 
{ T k\k=i °f users and any arbitrary chip waveform ip(t) satisfying the ICI free condition, the ESD of 
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the crosscorrelation matrix converges a.s. to the Marcenko-Pastur distribution with ratio index when 
M, K, N oo and K/N -> 0. ■ 

It is known that, in symbol-synchronous CDMA, the ASD of the crosscorrelation matrix at K, N — > oo 
and K/N — > /3 is the Marcenko-Pastur law [5, Proposition 2.1]. Thus, an equivalence result about 
symbol-synchronous and chip-synchronous CDMA in terms of the ASD's of crosscorrelation matrices 
can be established as follows. Under an unfaded channel 1 , the ASD of the crosscorrelation matrix in 
a chip-synchronous system converges, as M increases, to the ASD of the crosscorrelation matrix in a 
symbol-synchronous system with the same K/N ratio. 

To consider a more realistic scenario that the signal is subject to a fading channel, we define a quantity 
analogous to fJ,(R^) given in (24), expressed as 

n((AiR ct A) n )= lim {2M +ir 1 K- 1 E{tv((A^RipA (K Y)} , n>0. (26) 

M,K,N x I J 

K/N-+/3 

We will show below that (2M + l)' 1 K ~hr ((A^ R^A^) 71 ) converges to its limiting mean, i.e., 
»((A^R cs Ar). 

Lemma 3: Let denote the n-th moment of the random variable governing the asymptotic empirical 
distribution of the square magnitudes of received amplitudes {|^4fc(m)| 2 : — oo < m < oo,k = 
1, 2, ■ ■ ■ , K}. When M, K, N -> oo with K/N -► 0, the n-th moment of the ESD of A^ R ( c f } A^ 
converges to the nonrandom limit fi((A^ R cg A) n ), given by 

^R cs A r) = ±P^ £ n(n " 1) ^ (n " J \ +2) IT^H n>0, (27) 

j=l Cl+C2+ ... +Cj=n J{Cl,C 2 ,- ■ ■ ,Cj) r=1 

Cl>C 2 >--->Cj ->1 

where /(ci, C2, ■ ■ ■ , Cj) is defined in (62) of Appendix I. 

Proof: See Appendix IV. ■ 

We call fi( (A'R, s A)" ) the n-th AEM of the random matrix A^R Cfi A. The convergence of ESD of 
chip-synchronous CDMA in a fading channel is stated in the following theorem. 

Theorem 3: In a chip-synchronous system, if the moment sequence {7?( n )} holds for the Carleman's 
criterion, then for any relative delays {t^}^ 1 and arbitrary ICI free chip waveform ip(t), the ESD of 
{A^R^A^ : K = 1,2,-..} converges to a nonrandom limit whose n-th moment is equal to 
n((A^R cs A) n ) when M, K, N -> oo and K/N 0. 



'in an unfaded channel, the matrix A in (9) governing the amplitude of the received signal is a scaled identity matrix. Thus, 
it is the matrix R cs that determines the performance of the system. 
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Proof: By similar arguments as in [27], it can be shown that Y^=i ("P^ 2 "**) " = oo is a 
sufficient condition for Y,n=i Ki Aj RcA) 2n y 1/{2n) = oo, which implies {fi((A^ R cs A) n )} determines 
a unique distribution. Thus, with Lemma 3, this theorem follows directly from Theorem 1. ■ 

It is shown in [11] that the counterpart of R cs A) n ) in a symbol-synchronous system has the same 
expression as (27). Thus, in a fading channel, the ASD of the chip-synchronous system for large M is 
identical to that of a symbol-synchronous system, and the equivalence result of symbol-synchronous and 
chip-synchronous systems presented above for the unfaded channel can be generalized to the case of 
fading channel. Actually, the equivalence of the two systems can be discovered in an easier way. When 
all Tfc's are set to zero, a chip-synchronous system becomes symbol-synchronous. As Theorems 2 and 3 
hold for any realizations of relative delays, it is immediate to see the equivalence of chip-synchronous 
and symbol-synchronous systems. 

A related result has been demonstrated previously in [13]. By assuming the density of the relative 
delay distribution symmetric about NT c /2, it is shown in [13] that, as M — ► oo, a lower bound of the 
output SIR of the linear MMSE receiver for chip-synchronous CDMA attains that of the same receiver 
in a symbol-synchronous system. It is known that, given the linear model of a received signal, the 
MMSE achievable by a linear receiver, and hence the maximum output SIR, is dictated by the empirical 
distribution of the covariance matrix of the random channel matrix. It follows that our equivalence result 
on the ASD's of the crosscorrelation matrices of chip-synchronous and symbol-synchronous systems 
assures the equivalence of MMSE receiver output SIR in the two systems. Thus, the equivalence result 
we establish above holds in a more general sense, and neither an assumption about the distribution of 
relative delays nor a bound is employed. 

B. Chip-Asynchronous CDMA 

In computing the moments ^(iZ™ ), defined as (24) with R cg replaced by i? ca , the relative delays r^'s 
are regarded as either deterministic constants or random variables depending on the bandwidth of chip 
waveform ip(t). To be specific, it is known that, to satisfy the ICI free condition, the minimum bandwidth 
of ip(t) is 1/(2T C ) [20], which corresponds to the ideal Nyquist sine pulse. When ip(t) has a bandwidth 
of 1/(2T C ), Tfc's are treated as deterministic constants in the calculation of /j,(R™ a ); when the bandwidth 
of ip(t) is larger than 1/(2T C ), r^'s are taken as independent and identically distributed (i.i.d.) random 
variables whose density function possesses certain symmetry. The reason for this setting is due to the 
property of chip waveform presented below in Lemma 4. Thus, when the sine pulse is employed, r^'s 
are deterministic and the expectation taken in fi(R ca ) is w.r.t. random spreading sequences. If other chip 
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waveforms are used, resulting in a bandwidth lager than 1/(2T C ), the expectation in fi(R ca ) is w.r.t. both 
spreading sequences and users' relative delays. 

Lemma 4: Denote the Fourier transform of a real pulse ip(t) by 

/oo 
^(t)e- jQt dt. 
-oo 

Let 

/oo 
ip(t)if)(t - x)dt 
-oo 

be the autocorrelation function of ip{t). Define 

H^KToM^o 1 ) (28) 
= i?v(( n o - ni)T c + rj - ?7i)#v((ni - n 2 )T c + 771 - 772) • ■ ■ i^((n m _i - n )T c + ?] m _i - r/ ). 
We have the following results about ^({nj}™^ 1 ; {^j}™^ 1 ). 

1) For any n G Z and {Vj}]^ G R m , we have 

E ^(Kj^MS 1 ) (29) 

{ni,— .ra^-ijel-oo^]™- 1 

= I^WI^dO, m = l,2,-.., (30) 

ZTTlc J-tt/T c 

if the bandwidth of ip(t) is less than 1/(2T C ), i.e., = for ft > vr/T c . 

2) For any no G Z, 770 G R, and i.i.d. random variables {r/j}™!^ 1 satisfying E {cos(27rA;r/j/T c )} = 
for any nonzero integer k, we have 

E Ei^HCiWS 1 )} (3D 

{ni,- ,n m _i}6[— 00, oo]" 1 " 1 
1 /"°° 

= / l^(^)| 2m dO, m = l,2,..., (32) 

if the bandwidth of ip(t) is greater than 1/(2T C ). 
Proof: (Outline) This lemma can be proved by applying Parseval's theorem repeatedly for each 
summation variable m, n 2 , • • • , n m -\ in (29) and (31). Since the arguments of i?^(-)'s are cyclic, i.e., in 
the forms of (n - n{)T c + r/ - 771, (m -n 2 )T c + 771 -772, • • • , (n m _i - n )T c + rj m _i - 770, the complex 
exponentials due to Fourier transforms of time-shifted autocorrelation functions cancel each other. For 
the detail of the proof, see Appendix V. ■ 
Convergence of the ESD of {R^ : K = 1, 2, ■ ■ ■ } to a nonrandom limit when M, K, N ^ co and 
if/iV — ^ /3 is proved below. We define as the quantity given in (30) and (32), i.e., 



1 r°° 

— —r / \v(n)\ 2m dn, 



Wu' ] ~ —, / |*(.n)r"dfi. m- 1.2..... 
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Lemma 5: Consider a chip-asynchronous system whose quantity corresponding to the chip 

waveform exists for all m > 1. When the sine pulse is employed as the chip waveform, the relative 
delays r^'s are treated as deterministic; while if the bandwidth of the chip waveform is larger than 
1/(2T C ), then r^'s are viewed as i.i.d. random variables having E{cos(2irnTk/T c )} = for any nonzero 
integer n. For both short-code and long-code systems, when M,K,N — ► oo with K/N — > [3, /i(i2" a ) 
exists and is given by 

^^E^ 1 e Z~*-"-V 1 l) n <>> »^ ™ 

6i>6 2 >->6 n -i+i>l 

Proof: See Appendix II. ■ 
In the proof, when the bandwidth of ^>(f) is greater than 1/(2T C ), the formula of //(i2™) is obtained 
by means of the chip waveform property in part 2) of Lemma 4, which holds when distribution of r^'s 
has E {cos(27rriTfc/T c )} = for any nonzero integer n. A special case for this zero expectation is the 
uniform distribution in the interval [0, rT c ), r G N, which encompasses the symbol quasi-synchronous but 
chip-asynchronous system considered in [18]. Thus, as the equivalence in AEM leads to an equivalence 
in ASD, Lemma 5 provides with a proof for the conjecture proposed in [18] that the symbol quasi- 
synchronous but chip-asynchronous system has the same performance as a chip-asynchronous system. 

Theorem 4: Suppose that the chip waveform ip(t) has a finite bandwidth denoted by BW. If the 
sequence : n > 1} corresponding to ip(t) satisfies ^° =1 (w$ n) /2flw) = oo, then 

the ESD of {Rffl : K = 1, 2, ■ ■ ■ } converges a.s. to a nonrandom limit whose n-th moment is equal to 
fi(R^) when M, K, N -> oo and K/N -► /?. 



Proof: We rewrite W^ n) in (30) and (32) as 



2n 

d/, 



where / € S if ^(2vr//T c ) / 0. The measure of S is equal to 2BW. It is clear that |^(2vr//T c )/^| 2 
belongs to the space of integrable functions, and the set S is a measurable subset of real numbers with 
the Lebesgue measure. By a generalization of Holder's inequality [38], we have 

< (2BW) 1 ~ k / n k/H , l<k<n. 

Thus, we have the product of W^'s in (33) upper-bounded by 

b r /n 

(2BWy-"rr» \W^" } 

r=l r=l 



n—j+1 rt— jr'+l 



fl W? r) < n " n = {2BW) n -W^\ (34) 
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We use similar arguments as in [27] to show )} satisfies the Carleman's criterion. That is, we can 

bound fJ,(R^) by 

„(*») < wwr^ftiM' 1 E " < "- 1 >-° + 1 > 

3=1 V 7 b 1 +b 2 +-+b n 

6i>62>->6„_i+i>l 

= (2BW) n " 1 W^ ) ^ 



, i " 6l+62+ ... +bn _, +1=n /( & i'^,---,^ + i^ 



(3 \ J 1 1 /n\ / n 



< (2BWr- 1 wJ l) (l + ^) 2n . (35) 



So, 



f; Mier 1/(2n) > (1 + ^) " 2 £ « l) H ~ 1/(2n) = «, 

n=l ^ ' ra=l 

It follows tha the moment sequence {^(12™ )} determines a unique distribution. Besides, pursuing the 
same lines of the proof for Lemma 2 presented in Appendix III, we can show the n-th moment of the 
ESD of R { c p converges a.s. to n(R^ a ) when M, K, N -> oo and K/N -> /3. Thus, this theorem follows 
directly from Theorem 1. ■ 
We now consider the situation that the signal is subject to a frequency-flat fading channel. Define a 
quantity jJ,((A^ R ca A) n ) analogous to n((A^ R Cfi A) n ) of (26) by replacing R cs therein with il ca . We give 
the following theorem. 

Theorem 5: When M, K, N -> oo with K/N -► 0, the ESD of {A (x)t ^f : K = 1, 2, • • • } 

converges to a nonrandom limit whose n-th moment is given by 

j=l 6i+b 2 H \-b n -j+i=n 

bi>62>->6„_ j+ i>l 

- e c-) pk'iph -ai. 06) 

Cl+C2 + ... +C3= „/l l'°2, ,O n -j+lj/(Cl,C 2 , ,Cjj <=1 r=1 

Ci>C2>--->Cj>1 

if the sequences {7?( n ) : n > 1} and {>V^ n) : n > 1} satisfy J2n=i (^ {2n) >V^ n) /2Sw) 1/(2n) = oo. 

Proof: First, we prove the n-th AEM of A^R clL A, i.e., /i((v4^il ca v4) ra ), is given as (36). The proof 
follows the lines of Lemma 3's proof given in Appendix IV. It can be shown that fi((A' R^^A)" ) is 



October 6, 2008 



DRAFT 



17 



expressed as (cf. (94)) 

n{n-j)\(j - 1)! 



lim k-'y y y 

K/N^/3 



K,N,M^oo ^— ' ' ^— ' f(bl,bo,--- ,b„-i+l)f(ci,C2, 

j=l b 1 +b 2 +-+b n - j+1 =nc 1 +c 2 +-+c j =n J v 2 ' ' n V ^ ^' 



6l>6 2 >---6n- 3 '+l> 1 Ci>C 2 >-"C J >l 

j-1 n-j+1 j 

xY[(K-s)-N-* +1 11 ■ II^ (Cr) ' (37) 

s=0 t=l r=l 

which is equal to (36). 

Secondly, we would show Y^=\ {v^ 2n ^W^ n ^ /2B\vj ^ = oo is a sufficient condition that the 
sequence {n{{A^R c& A) n )} determines a unique distribution. By a generalization of Holder's inequality 
[38], we have V^ k) < (V {n) ) k/n for 1 < k < n. Consequently, the product of p( c -)'s in (36) is bounded 
as 

W. V {Cr) <[P (n) ) =V i - n) . (38) 

r=l 

Incorporating the inequality of (34), we can upper-bound ii({A' R^A)" ) by 
A* 



((AtJUA)") < (2BW)»- 1 wi B) ^ (B) E (J^) 

n(n- - 1)! 



x y y 

bl+b2+ .^ bn _ j+1=nCl+C2 ^ +C]=n f( b ^--- ,6„- j+ i)/( Cl ,c 2 ,--- , Cj ) 

fe 1 >6 2 >---b„„ :j+1 >l Ci>c 2 >---C3>l 

Proceeding in a similar way as the proof of Theorem 4, we are able to demonstrate that the condition 
Tin=i( vi2n) ^ n) / 2Bw y 1/(2n) = 00 is sufficient for J2n=i ^{{A^ R^A)^)' 1 /^ = 00, which 
gaurantees that {^((A^ R ca A) n )} determines a unique distribution. ■ 

We use the following corollary to establish the equivalence result of systems with three synchronism 
levels when M — ► 00 and the sine chip waveform is employed. 

Corollary 1: If M — ► 00 and the ideal Nyquist sine chip waveform 

^(i) = -^sincQQ (39) 

is used, the ESD of {Rffl ■ K = 1, 2, • • • } converges to the Marcenko-Pastur law with ratio index (5. 
Under the same premise, the ESD's of {A {K)] R<£) A {K) : K = 1, 2, ■ ■ ■ } and {^^i^fU^ : K = 
1, 2, • • • } converge to the same limiting distribution, provided that Y^=i (V^Tcj 1 ^ 2n ^ = 00. 
Proof: The Fourier transform of ip*(t) is 



= VT^rect 



2vr 
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where rect(x) = 1 for — 1/2 < x < 1/2 and equal to otherwise. By part 1) of Lemma 4, W^T = 1 
for all m G N. Due to (63), the formula of /t/(,R™ a ) in (33) is equal to (25), which is the n-th moment 
of the Marcenko-Pastur distribution. By the moment convergence theorem, the first part of this corollary 
follows. The proof of the second part is straightforward, where the equality of (96) is helpful. ■ 
It is demonstrated in [16] that, when the sine chip waveform is used and M — > oo, the asymptotic SIR 
at the linear MMSE detector output is the same for all of the three synchronism levels 2 . This equivalence 
result on output SIR can be seen as a direct consequence of the equivalence of ASD demonstrated by 
Theorem 2 and Corollary 1. It is shown in [39] that the linear MMSE receiver belongs to the family of 
linear receivers that can be arbitrarily well approximated by polynomials receivers 3 , i.e., in the form of 

f(A^RA) = a I + aiA^RA + ■■■ + a n (A^RA) n , 

with R standing for the crosscorrelation matrix in the system. In general, the accuracy of the approximation 
is in proportional to the order n of the polynomial. Both the coefficients a^'s and the receiver output 
SIR can be determined by the AEM of A^RA [6]-[9]. As AEM are equivalent in systems of three 
synchronism levels under the indicated circumstances, both the coefficients of the three polynomial 
receivers approximating linear MMSE receivers and their output SIR are identical. It is readily seen 
that the equivalence result is true not only for the linear MMSE receiver but also for all receivers in the 
family defined in [39], which proves the conjecture proposed in [16]. 

Up to now, the chip waveform is assumed to be ICI free for both chip-synchronous and chip- 
asynchronous systems. This ICI free condition requires that the chip waveform has a bandwidth no 
less than 1/(2T C ) [20]. Here we extend the equivalence results to the circumstance where the bandwidth 
of tp(t) is less than 1/(2T C ) so that zero ICI condition does not exist. At this moment, the crosscorrelation 
/9 cs (m, n; k, I) given in (7) is no longer correct. Instead, it has the same form as that of a chip-asynchronous 
system given in (13). The crosscorrelation in a symbol-synchronous system has the same expression 
as well by letting = t\ = 0. Setting ^'s in part 1) of Lemma 4 as the relative delays among 

users, it is shown by the lemma that J2{ ni ,- ,n m _ 1 }ehoo,oo]'™- 1 ^{{ n i}T=Q '■• i^}"^) 1 ) does not depend 
on realizations of relative delays. That is, this expression yields the same value in systems of three 



2 The equivalence results shown in [16] holds in a more general sense. That is, for any finite M, the output SIR of the MMSE 
detector in the chip-asynchronous system converges in mean-square sense to the SIR for the chip-synchronous system. 

'Although the result is presented in [39] for symbol-synchronous CDMA, the proof (Lemma 5 of [39]) can be extended to 
asynchronous systems in a straightforward manner. 



October 6, 2008 



DRAFT 



19 



synchronism levels. Tracing Appendix II for the proof of Lemma 5, we find out AEM formulas /u(i?" ) 
and fi((A^ R cs A) n ) have the same expressions as their counterparts in chip-asynchronous system, given by 
(33) and (36), respectively. Consequently, symbol-synchronous, chip-synchronous and chip-asynchronous 
systems have the same ASD when the chip waveform bandwidth is less than 1/(2T C ). Along with the 
equivalence result concerning the sine chip waveform in Corollary 1, the above discussion leads to the 
following corollary. 

Corollary 2: Suppose that M — ► oo and a chip waveform with bandwidth no greater than 1/(2T C ) 
is adopted. In either the unfaded or fading channel, systems with three levels of synchronism have the 
same ASD. ■ 

IV. More Results by Free Probability Theory 

In this section, we use free probability theory to obtain more results about the asymptotic convergence 
of eigenvalues of crosscorrelation matrices in asynchronous CDMA. Free probability is a discipline 
founded by Voiculescu [40] in 1980s that studies non-commutative random variables. Random matrices 
are non-commutative objects whose large-dimension asymptotes provide the major applications of the 
free probability theory. For convenience, the definition of asymptotic freeness of two random matrices 
by Voiculescu [41] is given below. 

Definition 1: [41] The Hermitian random matrices B and C are asymptotically free if, for all poly- 
nomials pj(-) and qj(-), 1 < j < n, such that /j,(pj(B)) = fi(qj(C)) = 0, we have 

f i(p 1 (B)q 1 (C)---p n (B)q n (C))=0. 

■ 

In this definition, the functional fi(-) is used. As we have shown in (24), p(-) is a limiting normalized 
expected trace of the matrix in the argument. Let B be sized by b x b, and we have a polynomial 

p( x ) = 2~27=o a i x% - Tnen 

n 

M(p(B))= lim b-^^EMB 1 )}. 

i=0 

Asymptotic freeness is related to the spectra of algebra of random matrices B and C when their sizes 
tend to infinity. In our context, the random matrices R cs and i? ca have column and row sizes equal to 
(2M + 1)K controlled by two parameters M and K. Since the asymptotes of ,R CS and ,R ca are studied 
when the size of observation window 2M + 1 is large, we let both M and K approach infinity. 
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In the following theorem, we show that R x , x G {cs, ca}, is asymptotically free with a diagonal random 
matrix D whose statistical description is detailed in the theorem. This asymptotic freeness property enable 
us to find the free cumulants of R x and AEM's of matrix sum R x + D and matrix product R X D. 

Theorem 6: Suppose that D_(m) = diag{c?i(m), d,2(m), ■ ■ ■ ,d,K(fn)} andZ) = diag{Z)(— M), D_(— M+ 
I))'"" iH(M)}, where dfc(m)'s are random variables having bounded moments, and dk{m\) and di{m,2) 
are independent if k ^ I. Also, R x , x G {cs,ca}, and D are independent. Then R Cfi and D are 
asymptotically free as M, K, N — > oo with i^T/iV — > /3. Moreover, if any of the following two conditions 
holds: 

1) The random variables c4(m)'s are non-negative, and exists for all m > 1, 

2) For any n G Z and {j?^}™^ 1 G M m , we have 

E l^tWSMCo^Ofl), m = l,2,... , 

ni,- ,n m _i£[- oo,co] 

then ^ ca and D are asymptotically free. 

Proof: See Appendix VI. ■ 
Before we proceed, some results of free probability theory about random matrices (see, for example, 
[42]) are summarized in the following theorem. 

Theorem 7: [42] Let B and C be asymptotically free random matrices. The n-th AEM of the sum 
B + C and product BC can be given by 

MOB + CD = E II (c\v\(B) + c lvl (C)) , (40) 

and 

M(Bcr) = EII c ivi( s ) II (41) 

Ve-ua U&KC{vj) 

where each summation is over all noncrossing partitions w of a totally ordered n-element set, V G w 
means V is a class of w, \V\ denotes the cardinality of V, Ck(B) is the /c-th free cumulant of B, and 
KC{w) is the Kreweras complementation map of w. Moreover, the relations between the asymptotic 
moment and free cumulant sequences are 

n{B n ) = EII C |V|( B )' ( 42 ) 

c n(B) = e n ^) n s \u\> ( 43 ) 

w V&vj U&KC(m) 

where 
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With the aid of Theorem 7, we consider free cumulants of R x and A*R X A for x G {cs,ca}. Rewrite (33) 
as 

Mny-E E n V" h " (n .7.) +2) n<^-'- w 

6!>6 2 >--->6 3 >l 

Let us interpret the summation variable j in (44) as the number of classes of a noncrossing partition w 
of an n-element ordered set, and b r is the size of the r-th class of w. From (42), it is readily seen that 
the n-th free cumulant of ,R ca is 



c n (R ca ) = W^^ n -\ 
Similarly, we obtain the n-th free cumulant of ,R CS as 

c n (R cs )=p n - 1 . 

Regarding the free cumulants of A*R CS A and A^R ca A, they are difficult to be identified directly from 
(42). Instead, we rewrite (43) in a more detailed way as 

n{n-j)\(j - 1)! 



= E E E 

j=l 6i+6 2 H \-bj=n ci+c 2 H hc n - j+1 =n 

bi>6 2 >->6j>l ci>c 2 >->c„_ 3+ i>l 



6 2 , • • • , b j )f(c 1 ,c 2 , ■■■ , c n _ j+1 ) 



j n—j+l 

xf[KB br ) ft S c t - (45) 

r=l t=l 

As AEM's are available for both A^R CS A and A^R ca A in (27) and (36), respectively, their free cumulants 
can be computed from (45). 

Let D be a (2M + 1)K x (2M + diagonal random matrix with the statistical properties stated 
in Theorem 6. Since D and R x , x G {cs, ca} are asymptotically free, (40) and (41) hold. Suppose that 
either the AEM or free cumulants of D are available. We have the n-th AEM of R ca + D and R C& D 
given as 

«t*.+m-± e , ^: l : ( :7: i| nK»-H(°^ <«> 



and 



6!>6 2 >->6,>l 



= E E E 



6i>62>->6j>1 ci>c 2 >->c„-i+i>l 



X 



i n-j+l 
r=l t=l 
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By setting D = AA^ , we have ^(D k ) = V^ k \ where is defined in Lemma 3. In this way, (47) 
becomes (36). 

The AEM ^((R Cfi + D) n ) and fi((R cs D) n ) can be obtained from (46) and (47), respectively, by setting 
all W^'s equal to one. In this way, /i((.R cs £)) n ) has a simpler form of 

n , .. x , . 1 n n-j+1 

„((^)»)=e^- j ' e '.'t +) n m^). 

j=l c 1+C2 +-+c„_ 3+1 =n^ Cl ' C2 ' ' C «-J+l^ r= i 

Ci>c 2 >--->c„_ J + i>l 

V. Connections with Known Results in Symbol-Synchronous CDMA 

We relate the results of this paper with those in [27], which find applications in symbol-synchronous 
CDMA. Consider a symbol-synchronous CDMA system. Define C = [c x c 2 • • • c K ] where c k is the 
N x 1 random spreading sequence vector of user k. Let S be an N x N symmetric random matrix 
independent of C with compactly supported asymptotic averaged empirical eigenvalue distribution. It is 
shown in [27] that the ra-th AEM of 



C T SC 

is given by 



(48) 



j=l b 1 +b 2 +-+b n - j+1 =n J y U 2 ' ' J + U r=l 

6 1 >b 2 >->6„- 3+ i>l 

We now establish the relationship of /i(-R™ ), /i(-R"a) ar >d n{{C T SC) n ). Denote the spreading sequence 
vector of user /c's m-th symbol as c k {m) = [c k mN ^ c ( mJV+1 ) . . . C (( m + 1 ) N - 1 )]T ^ an( j we jg^g 

C(m) = diag{c 1 (m),c 2 (m), • • • ,c K (m)}, (50) 
C = diag{C (-M),C(-M + 1), • • • , C(M)}. (51) 

Let A be a block matrix whose (to, n)-th block, — M < m,n < M, is denoted by A(m, n). Each 
A(m,n) is also a block matrix with the (/c,Z)-fh block, 1 < k, I < K, represented by A(m, n; k,l). 
The matrix A(m,n; k,l) is an iV x N matrix whose (p, q)-th entry, < p, q < N — 1, is equal to 
6((mN + p)T c + Tfc, (raiV + q)T c + 77). Then, the (to, n)-th block's (fc,/)-th element of R Cfi can be 
expressed as 

[Rcs]mn,kl = c k (m) T A(m, n; k, l)c t (n), (52) 
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and the crosscorrelation matrix R cs can be decomposed as 

R cg = C T AC. 

Similarly, we have 

[Rc a ]mn,ki = c k {m) T n(m, n; k, Z)cj(n) and R ca = C T (IC, (53) 

where matrix fl has the same structure as A with the (p, q)-th component of Q(m, n; k, I) equal to 

R^{{{m - n)N + (p - g))T c + T k - n). Rewrite n(R£) in (25) as 

j=l 6i+62+-+6n-i+i=n J V ^' ' " r=l 

bi>6 2 >->6n-j+i>l 

We find that fi((C T AC) n ) given m (54) and fi({C T QC) n ) given in (33) show remarkable similarity as 
fi((C T SC) n ) in (49). However, even though AEM's of matrices C T SC, C T AC and C T QC have the 
same form, they have distinct structures. As seen in (48), elements in the matrix C T SC are quadratic 
forms cj^Sci of a common matrix S. Whereas, in C T AC and C T ClC, the entries are expressed as (52) 
and (53), respectively, with the matrices A(m, n; k, I) and Q(m, n; k, I) varying for each component. 

In the following, another expression of R Cfi will be presented. Let u k (m), 1 < k < K and — M <m< 
M, be an A^-dimensional column vector whose y/N times scaled entries are i.i.d. random variables with 
zero-mean, unit variance, and bounded higher order moments. Besides, u k (m) and Ui(n) are independent 
when either k / I or m / n. Given a set of integers {-y k : 1 < k < K} G [0,N — 1] K , define the 
(2M + 2)iV-dimensional vector u k (m) 

u k {m) = [ 0,0, -j- ,0 ,u k (m) T , 0,0, -j- ,0 ] T . 

(M+m)N+"i k times (M—m+l)N— -y k times 

We also define a (2M + 2)iV x K matrix U_(m), given by 

U(m) = [ui(m),u 2 (m), ■ ■ ■ ,u K (m)], (55) 

and a (2M + 2)iV x (2M + matrix 

[/= [C/(-M),C/(-M + l),--- ,C/(M)]. (56) 

We have the following theorem, whose proof demonstrates that R cs can be expressed as U T U with certain 
choices of {7^} and « fe (m)'s. 

Theorem 8: For any {^ k } k=1 £ [0, N — 1] K , the ESD of the random matrix U T U converges a.s. to 
the Marcenko-Pastur distribution with ratio index [3 when M,K,N — > 00 and K/N — ► /3. Moreover, let 
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D be a diagonal random matrix as stated in Theorem 6. Then the n-th free cumulant of UDU is equal 

to fi(D n )p. 

Proof: Setting u k (m) as the spreading sequence vector of the m-th symbol of user k and 7^ := Tfc/T c 
in a chip-synchronous system, we have U_ T ( m )LL( n ) = [Rcs]mn and U T U = R cs - Thus, the first part of 
this theorem is a direct consequence of Theorem 2. 

For the second part, we have fi({UDU T ) n ) = j3 ■ fj,((R cs D) n ), written as 

t e "V, ,( H) 2 ' nx^. 

Ci>C 2 >--->C3>l 

By the moment-free cumulant formula of (42), the n-th free cumulant of UDU T is /j 1 (D n )j3. ■ 
Let us particularly use t/(m)* and E/* to denote the matrices U_(m) and {/, respectively, when = 
for all k's. Clearly, U^U* is a block diagonal matrix with each block C£(m)^C/(m)* a crosscorrelation 
matrix in symbol-synchronous CDMA. It has been derived in [11] that the n-th free cumulant of 
U.( rn )*D_(m)U_(m)J is n(D.(m) n )l3, which has the same form as its counterpart in chip-synchronous 
CDMA. 

VI. Spectral Efficiency and MMSE of Asynchronous CDMA 

In some applications of probability, it is frequent that the (infinite) moment sequence of an unknown 
distribution F is available, and these moments determine a unique distribution. Suppose that the final 
aim is to calculate the expected value of function g(X) of the random variable X whose distribution 
F is unknown. One of the most widely used techniques for evaluating E{g(X)} is based on the Gauss 
quadrature rule method [30], where 2Q + 1 moments {m n }^£ of distribution F are used to determine 
a Q-point quadrature rule {w q ,x q }® =1 such that 

/oo Q 
g(x)dF(x)^^w q g(x q ), 

and the approximation error becomes negligible when Q is large. However, this method often suffers 
from serious numerical problems due to finite precision of a computing instrument. Fortunately, by using 
the modified moments technique [43] which requires only regular moments {m n }, the algorithm becomes 
exceptionally stable especially when the density of the distribution F has a finite interval. In case that 
the interval is infinite, the algorithm does not completely remove the ill-conditioning [44, Section 4.5]. 
Some remedies can be found in the above reference. 

In this section, the Gauss quadrature rule method with modified moments technique is employed to 
compute the spectral efficiency and MMSE of asynchronous CDMA using AEM derived in previous 
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sections. Square -root raised cosine (SRRC) pulses with various roll-off factors a, denoted by SRRC-a, 
are adopted as chip waveforms. Since the employed method is numerically based and cannot refrain from 
computing errors, we are careful in drawing conclusions from the numerical results. Cares are taken to 
avoid making wrong claims caused by numerical errors. For example, that the spectral efficiency curve 
of system A is above the curve of system B may come from different amounts of numerical errors on 
spectral efficiency curves of the two systems. In the sequel, R ca , a (or R cs , a ) is used to represent the 
crosscorrelation matrix corresponding to the SRRC-a pulse. Given a random matrix M, we use Am to 
denote the limiting random variable governing eigenvalues of M when the matrix size tends to infinity. 

Assume the channel is unfaded and the per-symbol signal-to-noise ratio SNR is common to all users. 
We consider chip-asynchronous systems. The spectral efficiency of the optimum receiver is given as [2] 

C opt (a,/?,SNR) = JL E {log 2 (l + SNR- A^J} , (57) 

where the spectral efficiency is scaled by a factor (1 + a) -1 because the nonideal signaling scheme of 
SRRC-a pulse has each complex dimension occupies (1 + a) seconds x hertz. On the other hand, since 
the limiting distribution of the linear MMSE receiver output is Gaussian, the spectral efficiency of the 
receiver is asymptotically equal to the spectral efficiency of a single-user channel with signal-to-noise 
ratio equal to the output SIR of the MMSE receiver [2]. It is known that the MMSE receiver has the 
output SIR given as [5] 

I— -il- 1 
tr (I + SNRil caiQ ) -1, 

whose limit is lower-bounded by 

1 

1 + SNR • A^ 
and tr denotes the normalized trace. Thus, 

C mmse (a, 0, SNR) > --iLi g E J 1 I. (58) 

1 + a \ 1 + SNR • A^ J 

When a = 0, the equality holds 4 . To compare systems with chip waveforms of different roll-off factors, 

the spectral efficiency must be given as a function of the energy-per-bit relative to one-sided noise spectral 



4 The MMSE spectral efficiency can be obtained as C mmse (a,/3, SNR) = /31og 2 (l + SNR?7(SNR))/(1 + a), where t?(SNR) is 
the asymptotic multiuser efficiency of the linear MMSE receiver [45]. However, ?7(SNR) is not known to the author for nonzero 
a. 
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Fig. 1. Large system spectral efficiencies of a chip-asynchronous CDMA system for E b /No = 10 dB and unfaded channel. 



level Eb/No. It can be shown that a system achieving C opi (a, (3, SNR) has an energy per bit per noise 
level equal to [2] 



Eh 



/9SNR 



N (l + a)C°P t (a,/3,SNR)' 
and the same relation holds for the spectral efficiency of the MMSE receiver C mmse (a, (3, SNR) and 
E b /N . 

Fig. 1 shows the spectral efficiencies versus (3 in a chip-asynchronous system for the optimum and the 
linear MMSE receivers, where Eb/No is fixed as 10 dB and an unfaded channel is assumed. Spectral 
efficiencies corresponding to SRRC pulses with different roll-off factors are depicted, and curves in the 
figure are obtained from (57) and (58) using a 10-point quadrature rule. The black dots on the figure are 
obtained from the analytical results 

C°P\ a = 0,/3,SNR) = /31og 2 ^1 + SNR- ^(SNR,/^ + log 2 ^1 + SNR/3 - ^(SNR,/3)^ 

log 2 e 



4SNR 



:f(snr,/3) 



and 



C mmse (a = 0,/3,SNR) = /31og 2 1 + SNR - -^(SNR,/3) 
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when a = 0, where 

T{x,z) = (j + v^) 2 + 1 - ^(1- ^) 2 + l 

These results are derived in [2] for a symbol-synchronous system. However, by Corollary 1, they are 
applicable to a chip-asynchronous system with a = as well. It is seen that, when (3 is around 1, there 
is visible discrepancy between results of the analytical formula and the Gauss quadrature method on the 
spectral efficiency of the MMSE receiver. This is because the Marcenko-Pastur distribution, i.e., the ASD 
corresponding to a = 0, tends to be infinite-interval when (3 is close to 1, and the Gauss quadrature 
method is less accurate when the density function has an infinite interval. 

The discussions in the following two paragraphs apply to chip-asynchronous systems. For the optimum 
receiver, given any (3, the spectral efficiency corresponding to a = is obviously greater than that of 
a = 0.5 and then of a = 1. The spectral efficiency grows as (3 increases. When (3 is small, the ratios of 
spectral efficiencies of a = 0, 0.5 and 1 are roughly equal to the ratios of inverses of their bandwidths, i.e., 
ratios of (1+a) -1 , meaning the maximum bit rates that can be transmitted arbitrarily reliably are the same 
for various SRRC-q pulses 5 , although the consumed bandwidths are different. As (3 gradually increases, 
the ratios of spectral efficiencies (a = to 0.5 and to 1) become smaller and smaller, suggesting that, 
when a chip waveform with a larger excess bandwidth is chosen, the the maximum reliable transmission 
rate can be increased. 

For the linear MMSE receiver, the spectral efficiency is maximized by a certain (3 depending on a. 
When (3 is small, it is obvious that chip waveforms with smaller values of a have larger MMSE spectral 
efficiencies. Nonetheless, as (3 is greater than around 1.2, the favor of smaller a in spectral efficiencies 
disappears. For low (3, the linear MMSE receiver achieves near-optimum spectral efficiency. Otherwise, 
great gains in efficiency can be realized by nonlinear receivers. When (3 is small, the MMSE receiver 
with a = is superior to the optimum receiver having a = 0.5 in terms of spectral efficiency; so is the 
MMSE receiver with a = 0.5 to the optimum receiver having a = 1. Comparing curves of two receivers, 
we comment when more bandwidths are consumed due to the choice of an SRRC pulse with higher a, 
the return in channel capacity (maximum reliable data rate) is larger in the linear MMSE receiver than 
in the optimum receiver. For example, when (3 = 1, twice bandwidth of an SRRC pulse with a = 1 than 
with q = results in approximately twice reliable data rate of a = 1 than a = in the linear MMSE 



5 When a chip waveform with roll-off factor a is chosen, the maximum bit rates that can be transmitted arbitrarily reliably is 
equal to the spectral efficiency times (1 + ce)/T c . 
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Fig. 2. Large system optimum spectral efficiencies of symbol-synchronous, chip-synchronous, and chip-asynchronous CDMA 
systems when Eb/No = 10 dB and the channel is unfaded, where "ss", "cs" and "ca" refer to symbol-synchronous, chip- 
synchronous, and chip-asynchronous, respectively. 



receiver. However, for the optimum receiver, the ratio of the data rates between a = 1 and a = is 
around 1.5 for the same f3. Even for values of f3 in a practical system, the higher return of the MMSE 
receiver in capacity is still true. Another interesting observation is that the MMSE spectral efficiency 
curve of SRRC pulse with a = 1 is above the curve of the sine pulse in the region around (3 G [1.5, 2]. 
However, as the difference of the two curves is small, we should be careful in making comments. In 
that region, the curve of the sine pulse is exact (matching with the analytical result); the curve of SRRC 
with a = 1 is an approximation with two opposite forces counter-acting on each other. On the one hand, 
(58) is a lower bound so the curve underestimates the true spectral efficiency; on the other hand, the 
numerical method tends to be optimistic, which yields an overestimate. We surmise that the first factor 
dominates; adding the result that the curve of sine pulse is below the other, we conjecture: for some (5, 
the sine pulse is not optimal in terms of the MMSE spectral efficiency under the indicated environments. 

Fig. 2 shows the optimum spectral efficiencies as a function of (3 for symbol-synchronous, chip- 
synchronous and chip-asynchronous systems with various chip waveforms when E^/Nq = 10 dB and the 
channel is unfaded. The three curves marked by circles, squares and stars have appeared in Fig. 1 ; those 
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two marked by triangles (down and up) are obtained from (57), where Afl ca a in the equation is replaced 
with Afl cs a . In both symbol- and chip-synchronous systems, the spectral efficiency corresponding to a 
particular a* is equal to the spectral efficiency corresponding to a = divided by 1 + a*. This is because 
Afl cs a has the same distribution regardless of a. 

When a = 0, three systems have the same optimum spectral efficiency. This is a direct consequence 
of Corollary 2. Given any a, when (3 is low, the differences of the three systems in optimum spectral 
efficiency are negligible; as (3 increases, the chip-asynchronous system is superior to the other two for 
nonzero a, and the optimum spectral efficiency differences are proportional to [3. On the other hand, 
given any (3, as a increases, the chip-asynchronous system has a larger spectral efficiency than the other 
two, and the difference grows with a. Similar comments can be made from Fig. 1 for the MMSE spectral 
efficiencies of the three systems. We also observe, while choosing chip waveforms with larger bandwidths 
may result in the increase of channel capacity in a chip-asynchronous system, the statement is not true for 
symbol- and chip-synchronous systems. This can be interpreted as follows. It is the ASD that determines 
the performance measures of a system such as the channel capacity, the MMSE achievable by a linear 
receiver, and so on. Regretfully, the ASD of symbol- and chip-synchronous systems does not depend on 
the chosen chip waveform; hence the increase of bandwidth due to the replacement of a chip waveform 
merely decreases the spectral efficiency and does not help in boosting the capacity. 

Consider a channel subject to frequency-flat fading. The square magnitude of the received signal 
|^4fc("T-)| 2 is governed by SNR common to all users and a normalized random variable A^{m) having 
E{|Afc(m)| 2 } = 1. Thus, the amplitude matrix A has AA^ = SNRAA , where A has the same structure 
as the diagonal amplitude matrix A, and Ak(m) is located at the (k, fc)-th entry in the (m,m)-th block 
of A. The spectral efficiencies of the optimum receiver is given by (57) with A# replaced as A-rt D - T . 
Although we can also modify (58) to yield a lower bound for the MMSE spectral efficiency under fading; 
however, according to our experiments, the bound is loose. Fig. 3 compares optimum spectral efficiencies 
in a chip-asynchronous system with and without fading for a fixed E^/Nq equal to 10 dB. The fading 
channel is assumed to be Rayleigh. To generate curves of the fading channel, a 15-point quadrature rule 
is used. The black dots shown in the figure correspond to the analytical result obtained in [45] for a = 
and a Rayleigh fading channel. Perceptible discrepancy between analytical and numerical results appear 
in the region of (3 € [1,2]. We comment that fading in a chip-asynchronous system leads to a degradation 
in optimum spectral efficiency, which is consistent with a mathematical result demonstrated in [45]. 

In the presence of fading, the arithmetic mean over the users of the mean-square-error achieved by the 
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Fig. 3. Large system optimum spectral efficiencies of a chip-asynchronous system for unfaded and Rayleigh fading channels 
when E b /N = 10 dB. 



linear MMSE receiver is given as [5] 



For the unfaded case, £ (a, (3, SNR) is obtained by setting A in (59) as the identity matrix. Fig. 4 compares 
the MMSE achievable by a linear receiver for unfaded and Rayleigh fading channels when SNR = 20 dB. 
We use 10- and 15-point quadrature rules for unfaded and fading channels, respectively. Black dots on 
the figure correspond to the analytical result of 



for a = in the absence of fading [2]. According to our tests, the discrepancy between the analytical and 
numerical results grows with SNR. The difference becomes almost unnoticeable when SNR < 15 dB. For 
low (3, the MMSE in an unfaded channel is lower than in a Rayleigh fading channel. In the latter case, 
as [5 —> 0, the analytical result is equal to E{(1 + SNR • X)^ 1 } = 0.04079, where X has the exponential 
density e~ x , x > 0; the numerical result is equal to 0.03075 at (5 — > 0. 
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Fig. 4. Large system MMSE versus /3 of a chip-asynchronous system for unfaded and Rayleigh fading channels when 
SNR = 20 dB. 



We observe that, regardless of fading or not, the MMSE is inversely proportional to a in a chip- 
asynchronous system. This is consistent with the conclusion drawn previously that choosing chip wave- 
forms with larger excess bandwidths leads to a higher capacity. Nevertheless, for symbol- and chip- 
synchronous systems, regardless of a, the MMSE are unchanged and correspond to the curve of a = 0. 
Interestingly, we can see that fading decreases the MMSE in the region of high (3. The explanation is 
similar to that for fading increasing the spectral efficiency at high j3 made in [45]. That is, due to fading, 
a certain portion of interferers are low-powered; thus, the number of "effective" interferers seen by the 
receiver is reduced. This interference population control of fading compensates for its harmful effect on 
the desired user. It is also observed, as a increases, the receiver needs a larger j3 to have this phenomenon 
begin to operate, and this phenomenon is less obvious for larger a. 

VII. Conclusion 

In this paper, the ASD of crosscorrelation matrices in random spreading chip-synchronous and chip- 
asynchronous CDMA systems are investigated with a particular emphasis on the derivation of AEM. 
Noncrossing partition and the graphical representation of if -graph are the key tools in AEM computation. 
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We assume an infinite observation window width, known spreading sequences and relative delays to the 
receiver, and an arbitrary chip waveform. We consider both unfaded and frequency-flat fading channels. 
The spreading sequences are only assumed to be independent across users. For a particular user, we 
do not assume that the sequence is independent across symbols. Thus, results shown in this paper are 
applicable for both short-code and long-code systems. 

In the following, results of this paper are summarized. For chip-synchronous CDMA systems, the 
explicit expressions for AEM of the crosscorrelation matrix are given when the users' relative delays are 
deterministic constants. We show that AEM do not depend on the realizations of asynchronous delays 
and the shape of chip waveform, as long as the zero ICI condition holds. It is also shown that the 
AEM formulas are identical to those of symbol-synchronous CDMA. In an unfaded channel, as the 
AEM satisfy the Carleman's criterion and the a.s. convergence test, it is concluded that the ASD in a 
chip-synchronous system converges a.s. to Marcenko-Pastur law with ratio index 0. In the case of flat 
fading, the a.s. convergence of ESD to a nonrandom ASD is established provided that a constraint on 
the empirical moments of the fading coefficients is satisfied. 

For chip-asynchronous CDMA systems, the convergence of ESD to an ASD in a.s. sense is proved for 
general constraints on the chip waveform and, for a fading channel, the empirical moments of the signal 
received power. It is shown that, in contrast to chip-synchronous CDMA, AEM in a chip-asynchronous 
system are dependent on the shape of chip waveform. On the other hand, the relation of AEM and users' 
relative delays depends on the bandwidth of the chosen chip waveform. Specifically, it is mentioned that, 
for the zero ICI property to hold, the chip waveform has a bandwidth at least equal to 1/(2T C ), which 
corresponds to the sine pulse. When the bandwidth of the chip waveform is 1/(2T C ), AEM do not depend 
on the realizations of relative delays. On the contrary, if the bandwidth is wider than the threshold, AEM 
do depend on the asynchronous delays; nonetheless, different relative delays realizations may result in 
the same AEM. Suppose that relative delays are modeled as i.i.d. random variables Tfe's. Let G\ and 
C?2 be two distinct distributions of relative delays, and both of them possess the symmetry property of 
E{cos(27rriTfc/T c )} = for nonzero integer n. Then, for the same chip waveform, AEM's averaged over 
realizations of relative delays with distributions G\ and G2 are equal. The distribution symmetry condition 
given above encompasses the symbol quasi-synchronous but chip-asynchronous system considered in [18]. 
By moment convergence theorem, the equivalence of AEM leads to an equivalence of ASD provided 
that the uniqueness of limiting distribution is true. It follows that our result proves the conjecture given 
there that relative delays ranging uniformly within the chip duration and within the symbol duration 
yield the same performance. When the sine chip waveform is adopted, no matter fading or not, the AEM 
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of chip-asynchronous CDMA are shown to be equal to those of chip-synchronous CDMA and hence 
those of a symbol-synchronous system. This explains the equivalence result of [16] that the output SIR 
of the linear MMSE receiver converges to those of chip- and symbol-synchronous systems when M is 
large. Since every receiver in the family constructed in [39] can be arbitrarily well approximated by a 
polynomial receiver, and both the polynomial coefficients and the polynomial receiver's output SIR are 
determined by AEM, we can also prove the conjecture in [16] that the equivalence result in the output 
SIR of the three systems holds for all receivers in that family. We also study the situation that the chip 
waveform bandwidth is less than 1/(2T C ) such that zero ICI condition is lost. It is shown that, without 
the zero ICI property, the AEM formulas in symbol- and chip-synchronous systems bear the same forms 
as those in a chip-asynchronous system. Thus, when systems of three synchronism levels have the same 
parameters except for the delays of the users, their AEM's are all the same; consequently, these three 
systems have the same ASD. 

With the help of free probability theory, free cumulants of crosscorrelation matrices are also derived 
for both chip-synchronous and chip-asynchronous systems. It is also proved that the crosscorrelation 
matrix is asymptotically free with a random diagonal matrix having a general constraint. Based on the 
asymptotic freeness property, AEM's for the sum and the product of the crosscorrelation matrix and a 
random diagonal matrix are derived accordingly. 

Mathematical results obtained in this paper are connected to those that are widely used by researchers 
who apply random matrix theory to communication problems. 

At last, some application cases are provided. The Gauss quadrature rule method is adopted to depict 
the spectral efficiencies of the optimum and linear MMSE detectors and the MMSE achievable by a 
linear receiver in asynchronous CDMA. Performance in the measures of the spectral efficiency, channel 
capacity, and MMSE are compared for various chip waveforms, two types of receivers, and different 
asynchronism levels. 

Appendix I 
Noncrossing Partition 

The proofs of Lemmas 1 and 5 require results from noncrossing partition of set partition theory. Our 
treatment here for noncrossing partition is very brief; for more details, please consult [46]. 
Definition 2: (Noncrossing Partition [29,46]) Let S be a finite totally ordered set. 
1) We call w = {Bi, ■ ■ ■ , Bj} a partition of the set S if and only if Bi, ■ ■ ■ ,Bj are pairwise disjoint, 
non-empty subsets of S such that B\ U ■ ■ ■ U Bj = S. We call Bi, ■ ■ ■ ,Bj the classes of w. The 
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classes B\,--- ,Bj are ordered according to the minimum element in each block. That is, the 
minimum element in B^ is smaller than that of Bi if k < I. 

2) The collection of all partitions of S can be viewed as a partially ordered set (poset) in which 
the partitions are ordered by refinement: if w, a are two partitions of S, we have w < a if each 
block of w is contained in a block of a. For example, when S = {1,2,3,4,5,6,7}, we have 
{{1}, {2,5}, {3,4}, {6}, {7}} < {{1,3,4},{2,5},{6,7}}. 

3) A partition of the set S is called crossing if there exist pi < qi < P2 < 12 in S such that p\ and 
P2 belong to one class and q\ and q2 to another. If a partition is not crossing, then it is called 
noncrossing. ■ 

The set of all noncrossing partitions of S is denoted by NC(S). In the special case S = {1, ■ ■ ■ , n}, we 
denote this by NC(n). 

Definition 3: (Kreweras Complementation Map [29,46]) Consider elements 1,2, ■■■ , n and interlace 
them with 1, 2, • • • , n in the alternating way of 1, 1, 2, 2, • • • , n,n. Let w 6 NC(n). Then its Kreweras 
complementation map KC(vj) : NC(n) — > NC{n) G iVC({l, 2, • • • ,n}) is defined as the biggest 
element among those a G NC({1, 2, • • • , n}) such that the union zuUa of the two noncrossing partitions 
belongs to JVC({1,T, 2,2, ••• ,n,n}). U 
It can be shown that, if w contains j classes, then the number of classes in KC{w) is n — j + 1. 

A partition can be represented graphically. For example, Figs. 5(a) and 5(b) show two partitions of 
{ki,k2,--- ,kg}, where elements in the same class are joined successively by chords. A noncrossing 
partition is such that the chords intersect only at elements fci,--- , k n . For instance, Fig. 5(b) is a 
noncrossing partition, while Fig. 5(a) is not. In the following, we define a representation, called K- 
graph, for any partition of a totally ordered set. The K-graph defined below is similar to the VF-graph 
of [33] used to establish the convergence of moments of a Wigner matrix. We discover several pleasant 
properties of K-graphs that will be useful in proving the lemmas. They are enumerated right after the 
definition of K-graph. 

Definition 4: (K-graph) The i^T-graph corresponding to a partition w = {Bi,-- - ,Bj} of a totally 
ordered set {k\, ^2, • ■ ■ , k n } is denoted by a graph G = (V, £). The vertex set is V = {^1,^2, ■ ■ ■ , Vj}, 
and the edge set is £ = {e±, e2, • • • , e n }, where the edge e r connects vertices v s and vt if k r and k r+ \ 
are partitioned into classes B s and B t , respectively (with n + 1 := 1). ■ 

Remark: The fT-graph for a partition w of {fci, /S2, • • • , k n } can be interpreted in a more visually 
convenient way as follows. Let fcj's be arranged orderly (either clockwise or counter-clockwise) as vertices 
of an n-vertex cycle, and let edge e r connect vertices k r and k r+ \. The K-graph of w can be obtained by 



October 6, 2008 



DRAFT 



35 





2' k 7 , k t 



Fig. 5. (a) The partition {{fci, k±, ko}{k2, k^, kr, fcs}{&5}} of the totally ordered set {fei,fe, •■• , fcs}, (b) the partition 
{{fci}{A;2, /c3, fcr, fcg}{fc4, fcs}{fc6}}, (c) the if-graph for the partition represented in (a), and (d) the AT-graph for the partition 
represented in (b). 

merging vertices that are partitioned in the same class of w into one. When vertices are merged, edges 
originally incident on these vertices become incident on the merged one. Mergence of two adjacent 
vertices results in a self-loop cycle. ■ 
For example, Figs. 5(c) and 5(d) present the if -graphs for the partitions of Figs. 5(a) and 5(b), 



Properties of K- graphs: Given a partition w of a totally ordered n-element set K, and its corresponding 
if -graph G = (V,£). We have the following properties. 

1) There is a bijective correspondence between classes of w and vertices of G. 

2) G is connected. If and only if w is noncrossing, G is a concatenation of cycles with any two of 
them connected by at most one vertex. Moreover, if w is noncrossing and has j classes, then there 
are n — j '• + 1 cycles in G. For example, Fig. 5(d) is composed of 8 — 4 + 1 = 5 cycles. Any pair 
of these five cycles are connected by at most one of the two vertices labelled with k2,k^,kT, k$ 
and &4, &5. 

3) Consider a partition a of the ordered edge set £ = {ei, ■ ■ ■ , e n } of G by letting edges in the 
same cycle of G being in the same class. If w is noncrossing, then a is noncrossing as well. 



respectively. 
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Moreover, w and a are Kreweras complementation maps of each other. For example, Fig. 5(d) 
corresponds to a = {{ei, es}, {e2}, es, e$}, {e4},{e7}}, which is a noncrossing partition of 
{ei,e 2 ,--- ) e 8}- It is seen that {{1, 8}, {2}, {3, 5, 6}, {4}, {7}} is the Kreweras complementation 
map of {{1}, {2, 3, 7, 8}, {4, 5}, {6}}, and vice versa. ■ 

From properties 1) and 3), if w is noncrossing, then the corresponding i^-graph can represent both w 

and its Kreweras complementation map simultaneously. That is, w and KC (w) can be identified by the 

vertex set and edge set of the K-graph, respectively. 

Some results about noncrossing partition are in order. The number of noncrossing partitions that 

partition n elements into j classes is the Narayana number, given by 

1 f n\ ( n 



. , (60) 

Moreover, if the j classes have sizes ci , c 2 , • • • , Cj with c\ > c 2 > ■ ■ ■ > Cj > 1 (but not specifying 
which class gets which size), the number of noncrossing partitions is [29] 

n(n - 1) • • • (n - j + 2) 

77 \ , (61) 

/(Cl,C 2 , • • • ,Cj) 

where 

/(ci,c 2 ,-- - ,Cj) = Y[ n k- (62) 

k>l 

with nj; being the number of elements in (ci, c 2 , ■ ■ ■ , Cj) that are equal to k. It is clear to see 

n(n - 1) • • • (n - j + 2) = l/nW n \ 

Cl>C 2 >- - >Cj ->1 

The number of noncrossing partitions of an n-element set meeting conditions of 

i) w has j classes with sizes in non-ascending order of (ci, c 2 , ■ ■ • , Cj), and 

ii) the classes of KC{w) have sizes in non-ascending order of (pi, 6 2 , • • • , 6 n _j + i), 
is equal to [47,12] 



n(n-j)\(j - 1)! 



(64) 



f(h,b 2 , ■■■ , b n - j+ i)f(ci,c 2 , ■■■ ,Cj)' 

Appendix II 
Proofs of Lemma 1 and Lemma 5 

With the results of noncrossing partition in Appendix I, we now proceed to prove Lemma 1 and 
Lemma 5. Consider /u(i2™ ) in (24) and by replacing R Cfi in the equation with i? ca . They can be 
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rewritten as 

, x 6 {cs,ca}. (65) 



= lim 

K,iV —>oo 



lim (2M + l)- 1 E(tr((ilW) n )) 



For notational convenience, the superscript W of R[ K ^ will be omitted below when no ambiguity occurs. 
By (7) and (13), we have 

[Rcfi\m r m r +i,k r k r +i = Pes (jn r , ?7i r ^_i \ k r , fc r +l ) 
(m r +l)JV-l (m r+ i+l)AT-l 

E E C^C^SiPrTc + T^qr^n + Tk^), 

p r =m r .N q r+1 =m r+1 N 

and 

[-Rca]m r m r+ i,fc r fc r+ i = Pca(^V j "V+l j j ^V+l ) 
(m r +l)W-l (m r+ i+l)W-l 

E E ^ ) c^ ) ^((p r -g r+1 )T c + T fcp -7 lfcrfl ), 

p r =m r N q r+1 =m r+1 N 

for 1 < r < n with m n+ i := mi and fc n +i := fci. By expanding matrix multiplications of the term 
inside of square brackets of (65) can be expressed as 

lim (2M + l) _1 E{tr(i^)} 

M-+oo 

= lim (2M + V, 22 E {[^x]mim 2 ,S: 1 fe[Bx]m 2 m3,fei3 ' ' ' [Rx]m n m 1 ,k n k 1 } , (66) 

where K = {h,--- ,k n }, X = [1, K] x ■ ■ ■ x [1, K] = [l,K] n , M = {mi, - ,m n } and }> = 

V v ' 

n times 

[-M,M} n . Equation (66) is equal to 



i(»+>r'EE E - E « { C*'-*') (cS-'^) - - ■ («ft- > «a ,) ) } 

Kex MeyVieZi P n ez n 
x5(piT c + T kl ,q 2 T c + T k2 )5(p 2 T c + r k2 ,q 3 T c + T ks ) ■ ■ ■ 5(p n T c + T kn ,qiT c + r fel ), (67) 

when x = cs, and 

iim ( 2M + !>-'££ E ••■ E E { (* , 4? > ) (* , 4', ,) ) • • • } («) 

xE{i? v ,((pi - g 2 )^ c + r fcl - T k2 )R^,((p 2 - q 3 )T c + n 2 - T ks ) ■ ■ ■ R^{{p n - q\)T c + r kn - r fcl )}, 



when x = ca, where V r = {p r ,Qr} and Z r = [m r N, (m r + l)N — l] 2 for 1 < r < n. Owing to the 
tri-diagonal structure of R Cfi shown in (10), there are constraints \m r — m r+ i\ < 1 for 1 < r < n when 
(67) is considered. Moreover, as stated in the beginning of Section III-B, the relative delays r k s are 
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viewed as deterministic when the chip waveform is the ideal sine pulse and viewed as random when 
otherwise, the expectation in the second line of (68) can be discarded when the sine pulse is adopted. 

Computations of (67) and (68) can be executed by considering the equivalence patterns of elements in 
/C = {ki, &2, ■ ■ ■ , k n }. As equivalence relation and partition are essentially equivalent, the computations 
of (67) and (68) can be carried out with the aid of set partition theory, where K, is a totally ordered set 
with ordering k\ >- k2 >-■■■>- k n , and k r and k s are partitioned in the same class if and only if they 
take the same integer in [1,-K - ]. Note that the ordering k\ >- &2 >- ■ ■ ■ >- k n is just an arrangement of 
objects kis as an ordered set. It is different from the ordering of the values taken by summation variables 
h,k 2 ,--- ,k n in [1,K]. 

In the following, the summation Ylicex * n (^7) and (68) is decomposed into several ones using 
properties stated in Appendix I. Let ^CA'= [1, K} n such that each element xj = (xj(l), • • • , Xj(n)) 
in Xj corresponds to a j-class noncrossing partition of an n-element ordered set. We mean xj corresponds 
to a partition by that Xj(s) = Xj(t) if and only if the s- and £-th elements are partitioned in the same 
class in that partition. Moreover, let Xj(bi, h?,-, ■ ■ ■ ,b n -j+i), with b\ > 62 > ■ ■ ■ > 0n-j+i > 1> stand 
for the union of Xj's whose corresponding noncrossing partitions have Kreweras complementation maps 
with class sizes (61, 62, • • • , b n -j + i) (but not specifying which class gets which size). Since the Kreweras 
complementation map of a noncrossing partition is noncrossing as well, by (61), the number of elements 
in Xj(pi, 6 2 , • • • , b n - j+1 ) is given by 

62, • • • , 6n- J+ i) = ^~ h l) '" i i + l \ ■K(K-l)...(K-j + 1). 

/(Pi, »2, ■ • • > On-j+l) 

The above equation is interpreted as follows. The number of noncrossing partitions associated with 
Xj(pi, 62, ■ • ■ ,b n -j+i) is n(n — 1) ■ ■ ■ (j + 62, ■ ■ ■ ,6 n -j+i), and each of these noncrossing 

partitions has j classes, with each class specified by a distinct integer in [1, K]. Moreover, let X cro C X, 
where each element in X cro corresponds to a crossing partition of an n-element ordered set. With these 
settings, the summation Ylicex m (67) and (68) can be decomposed as 

n 

£ s E E E + E ■ m 

KeX j=l 6 1 +6 2 + ---+6„_ 3+1 =rt KeA^&i.&a,- ,b n - j+1 ) KdX cro 

6 1 >6 2 >--->6„_ 3+1 >l 

Now, we consider K-graphs corresponding to elements of Xj{b\, 62, ■ ■ ■ , On-j+i) and of X cro in (69). 
To embed the summation variables p r 's and qv's of (67) and (68) into a K-graph, in the n-vertex cycle 
composed of vertices k±, k2, ■ ■ ■ ,k n , two ends of the edge connecting k r and k r+ i are labelled with p r 
and q r +i, with the former and latter touching k r and k r+ i, respectively. We call these p r 's and q/s as 
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(a) 



(b) 



Fig. 6. 



The /("-graphs of (a) n-class partition {{ki} , {k 2 } , ■ ■ ■ , {fc„}}, and (b) (n 



l)-class partition 



{{fci}, {k 2 }, ■ ■ ■ , {k s , kt}, ■■■ , {k n }} with s < t. 

edge variables. Fig. 6(a) shows the graph with edge variables labelled. The edge connecting vertices k r 
and k r+ i stands for 



in (67) and (68), respectively, and the summands in the equations are yielded by multiplying terms 
associated with all of the edges together. Vertices fe r 's and edge variables p r 's, qy's in Fig. 6(a) are all 
summation variables of (67) and (68). Another set of summation variables m r 's are implicitly embedded 
in p r 's and qy's by that the ranges of p r and q r are both [m r N, (m r + 1)N — 1]. We would like to 
inspect under what equivalence relation of p r 's and g r 's would the contributions to (67) and (68) become 
nonvanishing in the large-system limit. 

A. Contributions of Noncrossing Partitions to [J>(R™ S ) and (j,(R™ a ) 

To compute the contribution of noncrossing partitions of K, = {k±, k<i, • • • , k n } to /i(-R™s) an d M(-^ca)> 
we replace Ylictx m (^7) an< ^ (^8) with the first term at the right-hand-side of (69) and then use (65). 
Given non-ascending ordered natural numbers (&i, &2> • " > &n-j"+i) sucn that 61 + 62 + • • • + & n -j+l = n - 
The contribution of each element of ^ (61, 62, • • • , ^n-j'+l) to (67) and (68) is to be evaluated. First, we 
consider j = n. There is only one if -graph, shown in Fig. 6(a). Since all of fci, &2, • • • , k n are distinct, 



4 ' S (Pr T c + T kr ,qr+lT c + T kr+1 ) 



and 




Rf((pr ~ q r +i)T c + r kr - r kr+1 ) 
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the expectations of spreading sequences in (67) and (68) are nonzero (equal to N~ n ) if and only if 
Pr = Qr for 1 < r < n. Note that, since c k Pr ^ and cjj?*^ are chips in the same transmitted symbol, i.e., 
m r N < p r , q r < (m r + 1)N — 1, the necessary and sufficient condition p r = q r stated above holds for 
both short-code and long-code systems. The contribution of each element of X n {n) to (67) becomes 

N~ n lim (2M + 1)- 1 V V V ■•• V 8{p 1 T c + T kl ,p 2 T c + r k2 ) 
Mey piez[ p 2 ez^, p„g.z; 

x 5(p 2 T c + r fc2 , p 3 T c + T k3 ) ■ ■ ■ 5{p n T c + r kn lPl T c + r kl ) , (70) 

where Z' r = [m r , (m r + 1)N — 1] and \m r — m r+ i\ < 1. The product of <5(-, -)'s is nonzero and equal to 
one if and only if all p r T c + r kr , 1 < r < n, are equal. Since we have pi G [miiV, (mi + l)iV — 1] for 
each mi € [— M, M], it is not difficult to see, in (70), the term behind N~ n is equal to N. Consequently, 
(70) is equal to N~ n ■ N = N~ n+1 . On the other hand, the contribution of each element of X n (n) to 
(68) is 

N~ n lim (2M + 1)- 1 E E E E{^((pi-p 2 )T c + T fcl -r fc2 ) 

M — ^oo — — — — 

X i?^((p 2 - P3)^ c + Tfc 2 - 75fc 3 ) • • • R^({p n ~ Pl)T c + Tfc„ - TfeJ} . (71) 

By Lemma 4, it is readily seen that, in (71), the term behind N 11 is equal to NW^ . Consequently, (71) 
is equal to N~ n • NW^ = N~ n+1 W^ . Note that the equivalence condition of p r = q r for 1 < r < n 
makes edge variables touching the same vertex in the K-graph, i.e., Fig. 6(a), to be equal. 

Next, we consider j = n — 1. Following the statement in the remark of Definition 4, we understand 
that any i^T-graph corresponding to <Y n _i(&i, can be obtained from the ivT-graph of X n (n) (denoted 
by G n ) by merging two vertices of G n into one. Suppose that vertices k s and kt of G n are merged and 
s < t, meaning that k s = k t in (67) and (68). Then, the original n-vertex cycle is decomposed into 
two cycles with numbers of edges equal to t — s and n — t + s as shown in Fig. 6(b). This K-graph 
corresponds to elements of A" n _i(max(t — s, n — t + s), min(i — s, n — t + s)). With a slight abuse of 
notational simplification, we let e = t — s and write A' n _i(max(t — s, n — t + s), min(t — s, n — t + s)) 
as X n -\(e,n — e) afterwards, although e may be smaller than n — e. 

We are going to demonstrate that, concerned with the four edge variables p s ,Qs,Pt,Qt of the merged 
vertex in Fig. 6(b), it is sufficient to consider p t = q s and p s = qt- Since fcj's are all distinct except for 
k s = kt, to yield a nonzero expectation of Y\r=i ^k^^k^ m (67) and (68), it is required that p r = q r 
for r e {1,2, •• • ,n}\ {s,t} and c^ k Pe \ c k 9s \ cj^ are in pairs, i.e., anyone of the following three 
conditions 
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Fig. 7. The directed flows represented by thick lines in (a), (b) and (c) denote the product of S functions corresponding to the 
equivalence patterns of i) q s — p s and q t = p t , ii) p t = p s + uN and q t = q s + uN and iii) q s — p t + uN and q t = p s — uN, 
respectively. 

i) q s = p s and q t = p t , 

ii) pt = p s + uN and q t = q s + uN with any integer u, 

iii) f/ s = ^ + uN and qt = Ps — uN with any integer it, 
for a short-code system, and 

iv) q s = p s and q t = Pt, 

v) p t = p s and q t = q s , 

vi) <? s = p t and q t = p s , 

for a long-code system. The reasons for same sign of two uN in condition ii) and opposite signs in 
condition iii) is because m s N < p s ,q s < {m s + 1)N — 1 and m t N < pt,qt < ("it + l)iV — 1. As 
conditions iv)-vi) ai - e special cases of i)-iii) with it = 0, we will use the latter to demonstrate our goal, 
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i.e., it is sufficient to consider only the equivalence relation of case vi). 

With conditions i)-iii) and p r = q r for r G {1, 2, ■ ■ ■ , n} \ {s, t}, the contribution of each element of 
X n -i(e,n — e) to (67) and (68) becomes 

A^+Mim^M + l)- 1 £ £ E E eR^*^} 

x5( Pl T c + r fel ,p 2 r c + T k2 )5(p 2 T c + r fc2 ,p 3 T c + r fe J • • • 5( Pn T c + r fcfi ,piT c + r kl )\ ks=kt , (72) 

and 

iv^+ 2 Hm (2M+1)- 1 v v E ••• E E R s) 4 9s) 4T t)c H 

xE{R^(( Pl -p 2 )T c + r fcl - T k2 )Rip({p 2 -P3)T C + r k2 - r fe J • • • 

Rf{(Pn ~ Pl)T c + T kn - T kl )} \ ks=kt , (73) 

respectively. A careful inspection reveals that the product of 5 functions in (72) becomes 

H S( Pl T c + T k ^, Pe T c + T k J, 2\ = {(1, 2), (2, 3), • • • , (n — 1, n), (n, 1)} (74) 

kt = k s 

for condition i), where we let q s := p s and q t := j>t, 

5(p s ^T c + T ks _^q s T c + T ks )5((q s + ^)T C + T ks ,p t -{T C + r fct _J 

X^fe + n^Te + Tfe^Pt+iT. + T^J 5{ Pl T c + T k ^p t T c + T kt ), 

(7,e)SX 2 

J 2 = {(1,2),-- - ,(5-2,s-l),(i + l,* + 2),--- ,(n,l)} 

U{(«,s + l),(« + M + 2),--- , (t - 2,t - 1)} (75) 
for condition ii), where we let pt := p s + uA r and g t := g s + uAT, and 

5(p s _iT c + T ks _„ (p t + uW) r c + Tk.)6(pt-iT c + n^ips - uN)T c + r fe J 



7 7"c + r k _ / , p e T c + r fce ) | J S(p v T c + r krt , p(T c + r k( ), 



k t =k s 

J 3 = {(1,2),--- ,{s -2,s -l),(t,t + !),■■■ ,(n,l)}, 

J 4 = {(a,a + 1),-" ,(t-2,t-l)} (76) 

for condition iii), where g s := p t + uN and g t := p s — uN. The product of functions in (73) 
can be obtained from the above equations by replacing each 5(a, b) with R 1 p(a — b). According to 
(74)-(76), conditions i)-iii) are graphically represented by Figs. 7(a)-(c), respectively. For instance, in 
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Fig. 7(a), except for the merged vertex, two edge variables touching the same vertex are the same. 
This is because p r = q r for r G {1, 2, • • • ,n}\ {s, t}, so we replace q r , r G {1, 2, • • • ,n}\ {s, t} in 
Fig. 6(b) with p r . Moreover, q s and qt in Fig. 6(b) are replaced with p s and p t , respectively, indicating 
that q s : = p s and q t := p t . In each of Figs. 7(a)-(c), the product of 5 and functions in (72) and (73), 
respectively, are depicted in the form of directed flow(s), indicated by thick lines traversing edges. The 
thick line passing through an edge with variables p 7 and p e (or q e ) represents 5{p 1 T c + r k ^,p e T c + r k J 
(or 5(p, / T c + T k _ t ,q e T c + r k J) for (72) and i?</,((p 7 - p e )T c + T ky - r fc J (or i?^((p 7 ~ <^) T c + T k-, ~ nj) 
for (73). For an edge with variables p 7 + uN and p e , it corresponds to 8{{p 1 + uN)T c + r ky , p e ^c + Tfc e ) 
for (72). Conditions i)-iii) are taken into account below. 

• Condition i): Consider the chip-synchronous case. The product of 5 functions in (74) is zero if 
p s / pt. Thus, we consider p s = pt, resulting in p s = q s = Pt = Qt, and the expectation in the first 
line of (72) is 0(N~ 2 ) because the fourth moment of is 0{N~ 2 ). It can be taken out from 
multi-dimensional summation, yielding 

0(N~ n ) lim (2M + 1)" 1 V V V V 6{p 1 T c + T kl ,p 2 T c + r k2 ) 

M — >oo L — * L — * L — * L — * 

M&yp 1 eZ' 1 p 2 &z 2 p n eZ' n 
xS(p 2 T c + T k2 ,p 3 T c + T k3 ) ■ ■ ■ 5(p n T c + T kn , Pl T c + r kl )\ k - kt . (77) 

It is readily seen that, in (77), the term behind 0{N~ n ) is equal to N. Thus, in this condition, the 
contribution of each element in <Y n _i(e,n — e) to (67) is 0(N~ n+1 ). 

Consider the case of chip-asynchronous. Assume p s / pt- We replace each 5(a, b) in (74) with 
R^(a — b) and plug the resulting product of R^, functions back into (73). The expectation of chips 
in (73) is iV~ 2 , and it can be taken out from the multi-sum, resulting in 

N~ n lim (2M + 1)- 1 E E E E{^((pi-p 2 )T c + r fel -r fc2 ) 
M&yp 1 eZ' 1 p 2 £Z 2 p n eZ' n 

X Ri>{{V2 ~ P3)T C + Tfc 2 - Tfc 3 ) • • • R^dPn - Pl)T c + T kn - T kl )} \ ks=kt . (78) 

By Lemma 4, the term behind N~ n is equal to NW^ for bandwidth of ip(t) either smaller or 
wider than 1/(2T C ). Thus, (78) is equal to N~ n+1 W^ . When p s = pt, the expectation of chips 
in (73) is 0(N~ 2 ), and (78) is revised by changing N~ n to 0(N~ n ) and imposing a constraint of 
p s = pt- Under these circumstances, the revised equation is equal to 0(N~ n+1 ). Thus, with both 
p s = pt and p s / pt, the contribution of each element of X n _i(e,n — e) to (68) is 0(N~ n+1 ). 

• Condition ii): Consider the chip-synchronous case. The product of 5 functions is given by (75). 
It can be checked that the product is zero if either u / or p s / q s not satisfied. Thus, we let 
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u = and p s = q s , which yields p s = q s = p t = qt- By similar arguments as in condition i), the 
contribution of each element of <Y n _i(e, n — e) to (67) is 0(N~ n+1 ). 

In the case of chip-asynchronous, we replace each 5(a,b) in (75) with R^(a — b) and plug the 
resulting product of R^, functions back into (73). For either p s = q s or p s / q s , the expectation 
of chips in (73) is 0(N~ 2 ), and we take it out from the multi-sum. Thus the leading term in (73) 
becomes 0(N~ n ). Tracing the proof of Lemma 4, we can find that the limiting normalized infinite 
sum of products of functions in (73), i.e., the term behind the leading term, is equal to 

N 



2ttT; 



— J e-^ 2uNT ^\^(n)\ 2n dn. 



When N is large, the integral is nonzero (equal to 2itT™ 1 W^) only if u = 0. Thus, the contribution 
of each element of AT n _i(e,n - e) to (68) is 0(N- n+1 ). 

Condition hi): In chip-synchronous case, the product of 5 functions is given by (76). This product 
is zero if u ^ 0. Thus, we consider u = and q s = p t / qt = p s , which implies m s = wit, and 
(72) becomes 

N~ n lim (2M + 1)- 1 V \{ 6( Pl T c + T^, Pe T c + T k J 

mi,— ,m s _i,mt,m t+ i,— ,m n (■y,e)e'I?, 
Pi,- ,p„-i,Pt,Pt+i,- ,Pn fc,=fc s 



X 



II *(p,r c + 7*„p f T c + Tfc 4 ). (79) 



m s ,m a+ i,— ,m t _i (jj,£)eZ 4 
p s ,p s+ i,- ,p t _i 
rrt 3 =m t 



Note that, for any particular mj, the constraint of m s = m< results in only iV choices for p s G 
[m t iV, (m t + l)iV - 1] in the second line of (79). It follows that (79) is equal to N~ n+2 . 
In the case of chip-asynchronous, (73) now becomes 

N~ n lim (2M + l) _1 Ei 



M^oo 



^2 Rip{( Ps -l - Pt - UN)T C + T ks _ 1 - T ks ) Yl R i>((P~f ~Pe)T c + T ki -TfcJ 



mi,— ,m,-i,mt,mt + i,- ,m„ ( 7j£ )gX 3 

Pl,— ,Ps-l,Pt,Pt+l,— ,Pn fe t =fc 3 

^ R^((p t -i -p s + uN)T c + r kt _ 1 -T fc J R^((p v - P()T C + r kv -r k( ) y 

m a ,m 3+1 ,— ,m t -i (n,C)£T4 ) 



(80) 



(«J,C)€ 

p s ,p s+ i,- ,p t _i 
m s =m t 



If the sine chip waveform is employed, the expectation of (80) can be removed. By similar arguments 
as in condition ii), the infinite sums in the second and third lines of (80) are zero if u ^ 0. By part 
1) of Lemma 4, when u = 0, (80) is equal to N~ n • Nwf ■ NW { £~ e) = N- n+2 wfw^' e) . 
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When the bandwidth of ip(t) is wider than 1/(2T C ), the calculation is more involved. We also let 
u = 0. Note that, in (80), both the product of functions in the second and third lines contain 
Tk s . We can rewrite the equation by 

N~ n lim (2M + 1)- 1 

E E II M(Pl~Pe)Tc + 

■ ,m 3 _i,m t ,m t+ i,--- ,m„ U 7i£ ) e x' 




mi,-- ,m 3 _i,m. t ,TO t+ i,--- ,m„ ^( 7)e ) fej , 3 

Pl,— ,Pa-l,Pt,Pt+l,— ,Pn fc t =fc 3 



Tfc 3 



E E | II ^((Pi - pO t c + r fc, - n ( 

+i,-,m t _i l(r?,C)eX^ 



(81) 



m s ,m s 

Ps,Ps + l,- ,Pt-l 

m 3 =m t 



where Tg = X3 U {(s — l,t)}, X4 = X4 U {(t — the first expectation at the second line is 

w.r.t. Tfc s , and the remaining two are w.r.t. {Tjt 4 }" =1 conditioned on Tk e . By part 2) of Lemma 4, the 

(e) 

two inner conditional expectations in the second and third lines of (81) are equal to NW^ and 
NW^~ e \ respectively. Thus, identical to the situation when tp(t) is the sine, the contribution of 
each element in X n -i{e,n - e) to (68) is equal to N- n+2 W { ^W^~ e) . 

To sum up, when j = n — 1, the contributions of each element in Xj(e, n — e) to (67) are 0(N~ n+1 ), 
0(N~ n+1 ), and N~ n+2 for conditions i)-iii), respectively. Regarding the contributions to (68), conditions 
i)-iii) have 0(N- n+1 ), 0(N- ,l+1 ), and N- n+2 W { fw^~ e) , respectively. Thus, as N -> 00, it is 
sufficient to consider condition iii) since the corresponding contribution has the highest order of N. Note 
that we require u = to get the result. 

We reach the conclusion that, when j = n — 1 with k s and k t partitioned in the same class, the 
equivalence conditions of edge variables {p r , q r } are that p r = q r for r G {1, 2, • • • ,n}\{s,t}, q s = pt, 
and p s = qt- Observing Fig. 6(b), we find that the equivalence relation can be stated as: within each of 
the two cycles in the K-graph, two edge variables touching the same vertex are equal. 

We consider j = n — 2 below. Two cases are possible. One is that k s , kt and k u , k v are respectively 
in the same class (l<s<u<v<t<n) and all others are singletons, whose corresponding K-graph 
is shown in Fig. 8(a). The other is that, except that k s , kt and k u (s < t < u) are partitioned in the same 
class, all others are singletons, whose K-graph is given in Fig. 8(b). In Fig. 8(a), the edge variables 
p s ,q s ,Pt, qt should be paired, and so do p u ,Qu,Pv,Qv By similar arguments demonstrated above for the 
case of j = n—1, we can see the equivalence relations of Pt = q s , Ps = qt, Pv = Qu and p u = q v yield a 
contribution to (67) and (68) with the highest order of N. On the other hand, in Fig. 8(b), the six edge 
variables p s ,qs,Pt,Qt,Pu,Qv should be in pairs, and the equivalence relations of p s = q t , pt = q u and 
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Fig. 8. The AT-graphs of (n — 2)-class noncrossing partitions of the totally ordered set K, = {ki, fa, • • ■ , k„}: (a) ks, kt and 
k u , k v (1 < s < u < v < t < n) are respectively in the same class and all other elements in /C are singletons, and (b) k s ,kt, k u 
are in the same class and all others are singletons. 



Pu = Qs yield a highest order of N in contribution. 

From the discussions about n — 2 < j < n, the following rule can be drawn. Let {Gj} denote the 
set of all if -graphs corresponding to j -class noncrossing partitions of K, = {kx, k<i, ■ ■ ■ ,k n }. From the 
noncrossing condition, it is immediate that, for any particular Gj in {Gj}, we can always find a member 
Gj+i in {Gj+i} such that Gj is obtained from Gj+i by merging two vertices in the same cycle of Gj + \. 
When the two vertices are merged, the cycle where these two vertices originally locate is torn into two. 
Within each of the two newly formed cycles, to yield a highest order of N in the contribution to (67) 
and (68), edge variables touching the same vertex should be set to equal. This observation leads to the 
following lemma. 

Lemma 6: Given non-ascending ordered natural numbers b\, hi, ■ ■ ■ , £>n-j+l sucn that bi + 62 + • • • + 
b n -j+i = n - Suppose that x E Xj(b\, 62, • • • , ^n-j+i)- Let G(x) denote the corresponding ET-graph of 
x, and c\(x) and C2(x) stand for the contributions of x to (67) and (68), respectively. To yield a highest 
order of N in c\(x) and C2(x), in every cycle of G(x), edge variables touching the same vertex should 
take the same value in [l,iV]. Moreover, we have 

n-j+1 

Cl (x) = N~i +1 + 0(N- j ) and c 2 (x) = N^ +1 W$ r) + 0(N~j). 

r=l 

Proof: The first part of the lemma can be proved by mathematical induction on j using the 
observation stated in the paragraph preceding this lemma. To prove the highest-order term in c±(x) 
is iV~ J+1 , we note that the expectation of n"=l c k c k m $7) and (68) is equal to N~ n when the 
equivalence relations of edge variables are satisfied. We take this iV~ n out from the infinite sum. The 
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limiting infinite sum of products of 5 functions in (67) can be decomposed into a product of smaller 
limiting infinite sum of products with each corresponding to a cycle in G(x) (see (79) as an example), 
and each decomposed limiting sum of products has a contribution of N to (67). Since the number of 
cycles in G(x) is n - j + 1, the term with the highest order of TV in ci(x) is N~ n n"=i +1 N = 
The proof that the highest-order term in C2(x) is N~i +1 Y\™~( +1 follows trivially from that for 

ci(x). m 
By Lemma 6, each element of Xj(pi, &2, ■ ■ ■ , b n -j + \) has the same contribution of iV _J ' +1 + 0(N~i) 
to (67). Thus, the total contribution of Xj(b\, 62, • • • , &„_j+i) is equal to 

#*i(&l,&2,-" ,K-i+l)-(N-' +1 + 0(N-')) 
= m + IlV - r) ■ + 0(*r'». 

}{oi,o 2 ,--- ,O n -j+l) 
Hence, by (65), the total contribution of noncrossing partitions of /C to fi(R^) is 

n(n - 1) • • • (j + 1) 

/(&!, &2, • • • , b n -j + i 



lim K^T V n(n - 1) ■ ■ ■ (j + 1) g _ _ +1 q 

jfjV^oo ^ . . ^ f bi.bo,--- ,6„_^i) liV 



where the equality holds by applying change of variables to (63). Similarly, the contribution of noncrossing 
partitions of /C to M(i2™ a ) is given by 

K ^- i f e / t^ ) ::f + 'V n' y - r »- jy - ,+1 "ff'< ) 

6i>6 2 >->6„-i+i>l 
n / i\ /• , i\ n-j+1 

- e^- 1 e ,t; + , n <>■ <», 

6i>6 2 >--->6n-j+i>l 

B. Contributions of Crossing Partitions to n(R™ s ) and fJ.(R r c l a ) 

Let G be a If -graph resulting from a crossing partition of K, into j classes. The graph G can be 
decomposed into at most n — j cycles. For example, Fig. 5(c) can be decomposed into at most five 
cycles. Thus, the contribution of G to (67) is 0(/V~- n • N n -i)Y\> r Zl(K - r) = O(l). By (65), the 
contribution of crossing partitions of K. to ^(R^) is 

n 

lim K- 1 #Xcro(3) ■ O(l) = 0, (84) 
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where #X cro (j) denote the number of j-class crossing partitions. By (82) and (84), we complete the 
proof of Lemma 1. Similarly, it can be shown that the contribution of crossing partitions of JC to /i(.R™ a ) 
vanishes in the large system limit. Thus, we have shown /i(-R™ a ) is given by (83), which completes the 
proof for Lemma 5. 

Appendix III 
Proof of Lemma 2 

By (19), the n-th moment of the ESD of when M — > oo can be represented by a normalized 
trace operator, given as 

r r ({R ( c P) n )= lim (2M + l)- 1 K- l tr({R ( c P) n ). 

M— >oo 

Define 

v K =Tr((RlPr)-KK)- 

By Borell-Cantelli lemma [48], if 

oo 

Prob(|^| > e) < oo, Ve > 0, 

K=\ 

then v k —>■ in a.s. sense. Using Markov inequality that 

Prob(|vx-| > e) = Pmb(v 2 K > e 2 ) < E{v 2 K }/e 2 , Ve > 0, 

we can show the a.s. convergence of vk to by proving ^2^=1 E { v \} < °°- I n tne following, we will 
show that 

E [tF((^f V) - »(R£ )] < oo (85) 

for all n G N, and the superscript of will be omitted for simplicity. In the proof, for simplicity, we 
suppose that the spreading sequences are independent from users to users, and for a particular user, the 
sequence is independent across symbols. That is, a long-code system is assumed. However, by following 
similar arguments of the proofs in Appendix II, it is straightforward to extend the proof here to a short- 
code system. 
We have 

E{[tfTO-M^ s )] 2 } 
= e{[5TO] 2 }-[m^ s )] 2 

= lim (2M + iy 2 K~ 2 y^Q(m u --- ,m n ,m n+1 , ■ ■ ■ ,m 2ri ; ■ • ■ , k n , k n+1 , ■ ■ ■ ,k 2n ), (86) 

M— >oo ' — * 
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where 

Q(mi,--- ,m n ,m n+1 , - ■■ ,m 2n ;ki,--- ,k n ,k n+1 ,--- ,k 2n ) (87) 

E{[-Rcs]mim2,fcife[^cs]m2m3,fofe ' ' ' [-^cs]m„mi ,k n ki 

x [-Rcs]m.„ + im„ +2 ,A;„ + iA;„ + 2 [-^cs]m.„ +2 m„ + 3,fc„+2fc„+3 ' ' ' [-^cs]m 2 „m„ + i,fc 2 „fc„ + i } 
— E{ [Rcs]m 1 m 2 ,k 1 k2 [Rcs]m2m 3 ,k 2 k 3 ' ' ' [Rcs]m rl m 1 ,k n k 1 } 

}, (88) 

and the summation is over all — M < mi, ■ ■ ■ , m 2n < M and 1 < k±, ■ ■ ■ , k 2n < K. Moreover, we have 



m; — m 



< 1 for {mj}" =1 , and so do {mj}fl n+1 



We consider two n-element noncrossing partitions. One is for {&!,••• ,fc ra }, and the other is for 
{fc ra+ i, • • • , k 2n ). Suppose that there are j and I classes in the former and latter partitions, respectively. 
Assume j classes of noncrossing partitions of {ki, • • • , k n } take distinct values {u±, • • • , Uj} in [1, K], 
and they have sizes (ai, ■ ■ ■ , aj), respectively. Similarly, {k n+ \, • • • , k 2n } take values {v\, • • • , v{\, which 
have sizes ■ ■ ■ ,bi), respectively. 

First, consider the case that {u\,--- ,Uj} and {vi,--- ,v{\ have no common element, i.e., all of 

(p) 

ui,--- ,Uj,v±,--- ,vi are distinct. Due to independence of spreading codes c\ s, the first term of 
(88) (expectation of a product of 2n elements) is equal to the second term. Thus, Q(mi,--- ,m n , 
m n+ i, • • • , m 2n ; fci, • ■ ■ , k n , k n+ i, • • • , k 2n ) is equal to zero, and (85) follows trivially. 

Secondly, consider the situation that {u\, • • • , Uj} and {v\, • • • ,vi} have only one element in common. 
Without loss of generality, we let u\ = v\ = w. In this case, (88) is equal to 



ai 



e [ n c^ct^ n & Mii) & iMii) i (89) 

a(l)=l 7 (1)=1 



e i n 

v a(l)=l 




Oi 

„(Pl,a(l))_(9l,o.(l)) I C I TT Jsi,7(l)) 



x ri E f ft & ,aw) <& i,aw) ] -II E ( II c& i,7(0) cS i,7<0) j x product of S functions, 

i=2 \a(i)=l / i=2 \ 7 «=1 / 

where, for given z and a(z), Pj )Ct (j) and gj )Ct (j) are edge variables touching the same vertex within a cycle 
of the i^-graph, and so do r i>7 ^ and Sj )7 (j) with given i and The summation in (89) is over all 
Pi,a(i)> Qi,a(i)> « = !,••• , j, = 1, • • • , ^ and r ij7(i) , s i)7(i) , i = 1, • • • , Z, 7(7) = 1, • • • , h L . Equation 
(89) is nonzero if and only if all the following conditions are met: 
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1) for each 2 < i < j, elements in {p i:1 , • • • ,p ij0i ,9i,i, ■ ■ • , 9;,aJ are in pairs, 

2) for each 2 < i < I, elements in {r^i, • • • , r^, s^i, • • • , Si >bt } are in pairs, 

3) elements of {pi,i,--- ,pi, ai , 9i,i, ■ ■ ■ , Qi, ai , ri,i, ■ ■ ■ ,r^ bl , s 1:1 , ■ ■ ■ ,s lybl } are in pair, and some of 
Pi,i, • • • ,Pi,oi, • • • ,91,0! are in pair with elements of r 1A ,- ■ • , r ljbl , Si,i, • • • ,s ljbl , 

where we say members of a set are in pairs, if each element of the set can find odd number of other 
elements that take the same value. To have the largest cardinality of summation variables mi, ■ ■ • , m2„, 
Pi,a(i)' s > 9i,a(i)' s > r i,7(«)' s and Sj )7 (i)'s in (89), the pairing constraints listed in conditions l)-3) above 
should be as least as possible, which yields 

4) for each 2 < i < j, {p ijl: ■ ■ ■ ,p i>a ., q iA , ■■■ , q ij0/ .} is paired by, for each 1 < a(i) < a u p i>a ^ is 
only paired with ?ija(i) , 

5) for each 2 < i < I, {r iyU • • • , r iM ,s i:1 , • • • , s ijbi } is paired by, for each 1 < < b h r ij7(i ) is only 
paired with s ij7(i ), 

6) there exists a unique (6>, z/) G [1, cti] x [1, denoted by (6q, uq), such that pi^ is paired with ri )I/0 , 
and q 1A> is paired with s ljVo (or pi j9o paired with s ljVo and gi^ paired with r 1>1/0 ), 

7) for {pi ; i, • • • ,pi,ai,9i,i, • • • ,9i, ai }, Pi,a(i) is on ly P^red with q lMl) when a(l) / O , 

8) for {ri,i,--- ,ri ibl ,si 5 i, • • • ,Si, 6l }, r lj7(1) is only paired with s 1)7(1) when 7(1) / fo- 
under these circumstances, the summand of (89) without the product of 5 functions is 0{N~ 2n ). On the 
other hand, for each mi and m n+ i, the product of 5 functions in (89) summed over dummy variables 
mi, - ■ ■ ,m n , m n+ i, • • • ,m 2n , Pi, a (i)' s > 9i,a(i)' s ' r i,y(i)' s and s »,7(j)' s is 

0(N n - j+1 ■ N n - l+1 ■ N- 1 ) = 0(N 2n - j - l+1 ), 

where N n ~i +1 and N n ~ l+1 are because {k\,--- ,k n } and {k n+ \,--- ,k2 n } form j- and /-class non- 
crossing partitions, respectively, and N _1 is because conditions 6)-8) causes the cardinality reduced by 
one 6 . Thus, (86) is qual to 

K- 2 ■ 0(N~ 2n ) • OiK'^ 1 ' 1 ) ■ 0(N 2n - j ~ l+1 ) = 0(K' 2 ), (90) 

where 0(if J+i_1 ) is because ui,--- ,Uj,v\,--- ,vi are all distinct except for u\ = v\. Consider the 
infinite sum over K of (85). It is finite when (90) is summed over K from 1 to 00. 



6 This is the same as conditions i) and ii) of Appendix II-A, where we compute the contribution of an element of X n -i(e, n — e) 
to (67). In these two conditions, edge variables touching the same vertex in a cycle of a /sT-graph do not always take the same 
value. 
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Next, we consider the situation that {ui, • • • ,Uj} and {v±, • • • ,v{\ have t elements in common, where 
2 < t < min(j, I). Without loss of generality, we assume m = Vi = Wi, 1 < i < t. In this case, (88) is 
equal to 

E|ri E ( ft *' a(i)) ^r w) ft ^: Mi)) ^: Mi)) ) pd 

[i=l \ a (i)=l 7 (i)=l / 

-fl E ( ft E ft ^: Mi)) ^: Mi)) 

i=l \ Q (i)=i / \ 7 (i)=l 

x ft E ( ft cS i,oC4>) cS 4,a(i)) ] II E ( 11 c^ w) c^ (l,) ) x product of<5 functions. 

i=t+l \a(i)=l / i=t+l \7(i)=l / 

The pairing constraints listed in 4)-8) is one of the conditions that yield (91) nonzero and have the 
largest cardinality of summation variables. Thus, the contribution of the sum of products of 5 functions 
is the same as that of when {u\, • • • , Uj} and {v±, • • • , vi} have only one common element. That is, it is 
equal to 0(N 2n ~^ l+1 ). It is not difficult to see, when {ui, ■ ■ ■ , Uj} and {v±, ■ ■ ■ ,vi} have t common 
elements, the contribution to (86) is 

K- 2 ■ 0{N- 2n ) ■ 0{K j+l - f ) • 0(N 2n - j - l+1 ) = OIK-*- 1 ). (92) 

When (92) is summed over K from 1 to oo, it is finite. 
Thus, we have proved 

oo 

^E{[fF( J R™)- M TO] 2 }<oo 

K=l 

for all natural numbers n, meaning that tr(il" s ) converges a.s. to ) when K — > oo. Moreover, it is 
seen ^t(i2" s ) given in (25) is equal to the n-th moment of the Marcenko-Pastur distribution. It has been 
proved in [26] that the moment sequence of the distribution satisfies the Carleman's criterion. 

Appendix IV 
Proof of Lemma 3 

Let P = AA^. Matrix P has the same structure as A. The fe-th diagonal entry in the m-th block 
diagonal component of P is Pk(m) = \ Ak(m)\ 2 . We have 

tr{(A^R cs A) n } = tr{(R cs PT}, 
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and 

\Z.\J -r- , ' ' '' ' 

O 

K/N^/3 



lim (2M + lr'K^EitvdR^Pr)} 



= Km (2M + ir 1 K- i y y y y y ... y 

K/N _ J =1 bi+b 2 +-+b n - j+1 =nKeX j (b 1 ,b 2 ,- K-j+i) MeyV^Zi V n eZ n 

bi>6 2 >-6n-j+i>l 

Ej^W) ( c £ 2) 4! 3) ) • ■ ■ (^"^l! ) } E{^ 1 (m 1 )^(m 2 ) ■ ■ ■ ^(m,,)} 

xS(piT c + T kl ,q 2 T c + T k2 )5(p 2 T c + T k2 ,q 3 T c + r fc3 ) • • • 5(p n T c + T kn ,q{F c + r fcl ), (93) 

with |mj — m-j+ij < 1. Consider the contribution of an element x in Xj(pi, b 2 , • • • , 6 n _ J+ i), which has a 
K-graph composed of n—j + 1 concatenated cycles with the number of edges b±, b 2 , ■ ■ ■ , b n -j + \ in a non- 
ascending order. Suppose that this i^-graph is yielded by a j-class noncrossing partition of K, with non- 
ascending class sizes (ci, c 2 , ■ ■ ■ , Cj). As shown in (64), the number of elements in Xj(b\,b 2 , • • • , 6 ra _j + i) 
satisfying these conditions is equal to 

n{n-j)\(j - 1)! 

/(&1,&2, • • • ,b n -j+l)f(ci,C2,--- ,Cj)' 

For x, the limit of the corresponding E{Pk 1 (mi)Pk 2 (m 2 ) ■ ■ ■ Pk n (m n )} in (93) becomes W^^V^, 
where note that if k s and k t of K are partitioned in the same class, m s is equal to m t . 7 Thus, for either 
short-code or long-code system, (93) can be expressed as 

n(n-j)\(j - 1)! 



lim K^y y y 



J=15l+fe2+ ... +fe _ +1=nCl+C2+ ... +Cj=n /( fe l' & 2,--- ,b n - j+1 )f( Cl ,C 2 ,... , Cj ) 

fe 1 >fe 2 >---fe„_ 3 + 1 >l Ci>C2>-C^>l 

x JJ(K - s) • A^' +1 • [J 7^), (94) 

s=0 r=l 

=i>~' e n( "»^ 2 ' ( ""^ 2) n^' (95) 

j=l Cl+C2 + ... +C3=n i^l^2, r=1 

Ci>C2>-">Cj>l 

where we make use of the equality 



fei+feH hb„-j+i=n 

bi>6 2 >-6n- J + l>l 



7 An illustration for the note is in case iii) of Appendix II-A, where we compute the contribution of an element of X n -i (e, n — e) 
to (67). 
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Appendix V 
Proof of Lemma 4 

We prove parts 1) and 2) of this lemma in Appendices V-A and V-B, respectively. 



A. Bandwidth oftp(t) Less Than 1/(2T C ) 

Define 6 1 = n Q T c + r]Q-'qi, 9 2 = n 2 T c + r] 2 -rji, R^,i{x) = R$(-x + B-y), and R\jj j2 (x) = R^{x-9 2 ). 
Then R^((no — n\)T c + t]q — rji) is the sample of R\p t i(x) at x = n\T c . The discrete-time signal R^,\(n) 
obtained from the continuous-time R^^x) by a sampling period of T c is denoted by the same notation, 
but the argument is an integer indicating the sample index. 

By Parseval's theorem, a partial sum w.r.t. n\ in (29) can be given as 



R^({no - ni)T c + r) - r]i)R^({m - n 2 )T c + rji - r) 2 ) 

n 1 =—oo 
oo 

n 1 =—oo 

= ±-J DTFT{^ i i(ni)}DFTF*{^ )2 (ni)}(L;, 
where DTFT{-} is the operator of discrete-time Fourier transform (DTFT) with 

oo 

DTFT{x(n)} = ^2 x(n)e~ 



(97) 



(98) 



As Rjf, : i(ni) is the sample of R^{—x + 9\) at time x = n\T c , and the Fourier transform of R$(x) = 
ip(x) * ip(—x) is |^(r2)| 2 , we have 



1 00 - / 

DFTF{ J R v ,, 1 (n 1 )} = - £ e"^* *( 

k=—oo 



' u — 2irk 



where uj = VLT C . Consequently, (98) is equal to 



2vrT 2 5^ / 



oj — 2irk 



uj — 2irl 



dco, 



(99) 



where, since ^(fi) is bandlimited to tt/T c , only k = I = has nonzero integral. Thus, (97) is equal to 



/e - 1 T c 



{(n Q -n 2 )T c +ri -i)2) 



du. 



(100) 
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We consider the summation w.r.t. 712 in (29). Define #3 = n^T c + 773 — 772 and .6^3(2) = R^(x — #3) 
We have 

4 \ 

duj R^{{n 2 - n 3 )T c + rj 2 - 7/3) 



E - 

n 2 =— 00 \ 



D -i^((no-n2)T c +J7o-r? 2 ) 



-j^r{n Ta+ria-V2) 



4 ( E e^^, 3 (n 2 )^dcu, 



(101) 



where the summation inside the brackets of (101) is the complex conjugate of the DTFT of ^,3(712), 
given by 



T 



k=— oo 



cj — 27r/c 



(102) 



Plugging (102) back to (101), we can see that the integral is nonzero only when k = 0, which results in 



1 



duo. 



(103) 



2vrT c 3 J_ v 

In consequence, when the summations w.r.t. n\ and n 2 are taken into account, the result is given in (103). 
Continuing this process, we obtain the final result as 



1 r 

2vrT c m J_„ 

which is equal to (30) by setting Q = u/T c . 



2m 



duj, 



B. Bandwidth of ip(t) Greater Than 1/(2T C ) 

Suppose that the bandwidth of ip(t) is greater than a/(2T c ) but less than (a + 1)/(2T C ) for a G N. 
Using the equality E UtV {g(u, v)} = E v {E u {g(u, v)\v}} with E u {-|-} denoting the conditional expectation 

w.r.t. u, we can see that 

00 

E E r/!{^(( n o ~ n i) T c + f]o ~ m)Rip(,( n i - n 2 )T c + 771 - 772)1772} (104) 

ni=— 00 

is nested in the multi-dimensional summation of (31), where note that 7/0 is deterministic. By employing 
the same procedures of getting (99), it is immediate (104) becomes 



I 00 PIT . 

27r T2 E / Er li{' 
171 ±c k,l=-oc J -* 







2 



















do;, 



whose imaginary part is definitely zero. Since {cos(27r/c77i/T c )} = for any nonzero integer k, it is 
readily seen that we only need to consider k = I in the above equation, which is given by 



1 r«/2i . w 



_ J -u ! ^r*(( no _ n2 ) Tc+J?0 _ J?2 ) 



fc=- fa/21 



a; — 27rA; 



dw. 



(105) 
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Consider one more layer of summations in (31), i.e., w.r.t. ri2, which yields 



1 fa/21 
^f2 XT 



2vrT 2 



fc=- fa/21 



\ 4 / oo 

' \?12= — OO 



dw, 
(106) 



where the term inside the brackets can be written as 

( oo 

E* e -^("o r '^-*) £ ^i^, 3 (n 2 ) „, 

I. ri 2 =— oo 

1 oo 

= J_ ^ E?? | e -i i ^K'j r <=+ 1 ?o-r ?2 ) e i^i(r i 3T c + J ?3- J 72) 
t=— oo 

with the expectation in (107) being nonzero only when I = k. Thus, the summations w.r.t. m and n 2 of 
(31), i.e., (106), become 



— 2"7r/ 



(107) 



1 r«/ai 

c fc=- fa/21 J 71 



.jS!_^^ no _ n3 ^ Tc+rlo _ rl3 ^ 



lo — 2nk 



dto. 



Continuing this process, we obtain the final result as 



1 ^ r f uJ - 27Tfc \ 



2m 



duo 



1 r (2\a/2\+l)n 



k=- fa/21 

which is equal to (32) by changing variable from w to Q = lo/T c 



2vrT c m y_(2[ Q / 2 i+i) 7r 



to 



2m 



d(J, 



Appendix VI 
Proof of Theorem 6 

The following two lemmas are helpful to prove Theorem 6. 

Lemma 7: Suppose that w is a noncrossing partition of a finite totally ordered set S, where every 
class in w has at least two elements. Then, there exist some classes in w that contain adjacent elements 
of S, where the adjacency is cyclic ordering, i.e., the first and last elements of S are adjacent. 

Proof: Denote the i^-graph corresponding to w by G{w). From the properties of a i^-graph given 
in Appendix II, there is a bijective correspondence between the class set of w and the vertex set of G{w). 
We have the following three observations. First, the size of the r-th class of w is equal to d(v r )/2, where 
v r is the vertex in G(w) that corresponds to the r-th class, and d(v r ) is the degree 8 of v r . Secondly, 
whenever two adjacent elements of S locate in the same class of w, there is a self-loop in G{w). Thirdly, 



The degree of a vertex is the number of edges that connect to that vertex. The singly vertex in a self-loop has the degree 
equal to two. 
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for all noncrossing partitions w of S, we cannot find any G(zu) that contains no self-loops and whose 
every vertex has the degree equal to or greater than four. Based on the first two observations, the third 
one can be interpreted as the statement of the lemma. Thus, we have completed the proof. ■ 
Lemma 8: Let D be the diagonal random matrix described in Theorem 6. We have, for x 6 {cs,ca}, 

lini (2M + iy 1 K- 1 £ E E {^x 1 ]m 1 m 2)fel foP Sl ]m 2 m 2 ,i fc2fe2 

[-^x ]m2m 3 ,k 2 k 3 ]m 3 m 3 ,k3k 3 ' ' ' [-^x \m n m 1 ,k n k 1 \D ]mim 1 ,k 1 ki} i (108) 

= p{B?) lirn (2M + 1)- 1 K- 1 £ £ 

K/N-+/3 K.\{k t+1 } M\{m t +i} 

E{(d k2 (m2)) Sl ■ ■ ■ (dfc t _i("»*-i)) at " !, (rffc t ("H)) at " 1+it (dfc, +2 ("»t+2)) at+1 • • • (4>i)) 5 "} 

m 1 m 2 ,k 1 k 2 ' ' ' [-R x * 1 ]m t -im t ,kt-ik t [Rx +1 \m t m t+ 2,ktkt+ 2 ' }, (109) 

where r^'s and Sj's are non-negative integers, notations /C, A', .M and ^ have the same definitions as in 
(67) and (68). When D is set as the identity matrix and the constraint k t = kt+i in (108) is replaced 
with kt = K for any t < u < n, (108) is equal to /i(^ t+ - +r "- 1 )^( J R^ + - +rt - 1+r " + - +r "). 
Proof: We expand the multi-sum of matrix product in (108) as 

e{k 2i1 (d U3A Kor • • • (d ulA my-} 

l<Mj,iU) <K,-M<Vj,,u) <M 
l<j<n,l<l(j)<rj 

U t ,l=U t +l 

xE{[R x ] 

fl,l"l,2,«l,l«l,2 Vl,2Ul,3>Wl,2Wl,3 Ul,ri«2,l,Ul,rilt2,l 

V v ' 

[iix \m 1 m 2 ,k 1 k 2 

}, (HO) 

v 1 1 1 1 — — ' ' '■ ' ' 

[R\ ]m n m 1 ,k n k 1 

where we let Uj t i := kj and Vj t \ := rrtj for 1 < j < n. To compute (110), we consider noncrossing 
partitions of the ordered set {uj^ : 1 < j < n, 1 < /(j) < r^} with a constraint that u tj i and ut+i t i are 
in the same class. The ordering of the set is clear from the equation. The .KT-graphs corresponding to these 
noncrossing partitions can be generated by vertex mergence of a -KT-graph, called original K-graph, similar 
to the one in Fig. 6(b). In specific, consider a cycle with r± + ■ • ■ + r n vertices. In a counter-clockwise 
direction, the vertices on this cycle are U = • • • , ui,n , ^2,1, ■ ■ ■ , ^2,r 25 " " " > u Hi i, • • • , u n ,r n }- The 

original K-graph is obtained by merging vertices u t: i and 1^+1,1 into one vertex. Let us call the cycles 
at the left- and right-hand-side as L-cycle and i?-cycle, respectively. 

By the property of noncrossing, it is sufficient to consider only K-graphs obtained by executing vertex 
mergence individually on L- and i?-cycles. That is, noncrossing partitions are performed respectively 
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on the two ordered sets {1/4,1, • • • , ixt, rt } and U \ {ut >2 , • • • , ixt, r , , ixt+1,1}. Note that the two sets have a 
common element ut : i. The .fT-graphs yielded in this way give non-vanishing contributions to (109) in 
the large-system regime. It follows that (110) can be written as 

£ E{(d U21 (v 2 ,i)) Sl ■ ■ ■ (d Ut _ 1A (v t - lyl )) St - 2 (d Ut l (v t ,i)) St - 1+St (d Ut+21 (v t+2 ,i)) St+1 

U\{u tt2 ,--- ,M t ,r t ,Mt+l,l} 
V\{t) t ,2,--- ,V t ,r t ,V t +l,l} 



(d^M^r-miR: 1 ] fi,if2,i,«i,i«2,i ["x Jft-i,ift,i,«t-i,iu t ,i L-"": 



X }Vt,lV t + 2,l,Ut,lU t + 2,l 



Vt,lV t ,2,Ut,l u t,2 L-"->tJft,r t Ut,l,Wt,r t Ut,l }, (HI) 



Vt,2,--- ,V t ,r t 



where V = {vi,i,--- , v 1<ri , v 2 ,i,- ■ ■ ,v 2 , r2 ,--- ,v n ,w- ,v n ,r n }, each [Rx]v i , lVi+1A ,u i , 1 u i+1A in the second 
and third lines is expanded as the product [■fl x ]u iil D ii2 ,u iil u ii2 [-ftx]u ii 2i; i ,3,rt i ,2u i ,3 '"' Ui.rj «i+l,l,«i,ril*i+l,l ' 
The integers in [1, K] chosen by elements in U \ {ix*,i, ixt,2, ■ ■ ■ , ut, rt ,ut+i,i} are all distinct from those 
chosen by elements of ut :2 , ■ ■ ■ ,ut^ t - That is, besides u^i, elements in sets {1*4,1, • • • ,ut,r t } and U \ 
{ut, 2 ,-" ,ut : r t ,ut+i : i} choose no common integers. So do sets V \ {«t,i, vt,2, ■ • • ,vt, rt ,vt+i,i} and 
^t,2, • • • )^t,r t - Note that, although the second and third expectations of (111) are both concerned with 
common summation variables u t: i and v^i, the random variables indexed by u t) \ and v^i can still be 
put in two different expectations. For details, see the discussion of condition iii) in Appendix II-A. Since 
the limit of the multi-sum in the third line of (111) is equal to n{R r ^) for any integers i*t,i G [1, K] and 
vt,i G [— M, M\, it can be factored out to the head of the equation, and we obtain (109). 

When D is set as the identity matrix and the constraint k t = k t +\ is replaced with k t = k u , (109) can be 
revised accordingly. We can see the revised equation is equal to /x(12£' H hru ~ 1 )/j,(R^ 1 ^ hr * _1+r, " H hr "). 

■ 

Proof: [Theorem 6] 
Suppose that, for 1 < j < n, polynomials 



Pj( x ) = a J>X' and = S b i> s * 

Tj>0 s 3 >0 

give 

rj>0 Sj>0 

We have 

(i( Pl (R x )q 1 (D)--- Pn (R x )q n (D)) 



S-i 

X 3 



V ai, ri 6i, Sl ---a n , rn 6„, Sn lim (2M + l)" 1 K^E^K^D^ ■ ■ ■ R^D S -)}, (112) 

z — * M,N,K^oo 



+00 

ri '"' ,r " K/N->[3 
si,-,s„ ' 
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where 



E{tr(il^£) Sl ■■■R r x "D s ")} 

^ ] E{ [i? x 1 ] mirrl2j fc 1 fc 2 [X) 1 ] m2m2) fc 2 ); 2 • • • [-R x "]m.„mi,fc„fci [-O "]mimi,fcifa} (113) 



mi,--- ,m„ 
fci,--- ,fc„ 



mi,-- ,m„ 

fcl ,'" ,^T1 

XE{ [^ x 1 ] miTra2j fe 1 fe 2 [^ x 2 ] m2rrl3) fc 2 fc 3 • • • [-R x n ]m n mi,fe„A;i} (114) 

E{(d U2il (u^))* (e^ ( W3 ,i)) S2 • • • 

1<«3,((J) < K - M < V l,l(j) < M 

l<j<n,l<l(j)<r 3 

><E{[R x ] Vl lVl 2tUl lUl 2 [R x ] Vl ^ Vl ^ 3jUl2Ul3 ■ ■ ■ [-Rx]j;i lt . 1 ll2,l,Wl,riW2,l 
V * ' 

[Ry. ]rn 1 m 2 ,k 1 k 2 

" [-Rx]ll„,iU„,2,U„,iU„,2 [■^x]«n,2fi.,3,Wn,2Un,3 ' ' ' ["^x] t!„ , r „ II i , 1 ,«„,,. lii, 1 } ■ (H5) 
V . ' 

[Rx n ]m n m 1 ,k n k 1 

In (115), we use Vj : \ := rrij and u^i := kj for 1 < j < n. Our goal is to show /j,(pi(R x )qi(D) ■ ■ ■ 
Pn {R x )q n {D)) = 0. 

To compute (115), noncrossing partitions of the ordered set {uj^ : 1 < j < n, 1 < Z(j) < rj} (the 
ordering is as shown in (115)) are considered. The partitioning can be decomposed into two stages. In 
the first stage, we perform noncrossing partitions on the ordered set {uj : \ : 1 < j < n}\ in the second 
stage, elements in {ujUj\ : 1 < j < n, 1 < < rj} \ {uj : \ : 1 < j < n} are partitioned into classes of 
the noncrossing partitions performed in the first stage according to the noncrossing condition. We divide 
the noncrossing partitions performed in the first stage into two groups as follows. 

Group 1: At least one of the classes contain only one element. 

Group 2: Every class contains at least two elements. 

For Group 1, without loss of generality, we suppose that u\ : i is a singleton. Then, the expectation of 
d Uj i(^j,i)'s in (115) can be written as 

ElK.J^^Oj^lElK^K!)) 5 ^^,!^,!))^---^,!^,!)) 5 "- 1 }, (H6) 
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since dk(mi) and di(m,2) are independent if k ^ I. Then, (115) can be written as 

£ E{(d UlA ( Vl ,i)) Sn } x 

Ul,l,Vl,l 

2 E{« il Ki)) Sl K 3 , I Ki)r-(4„, 1 (M)) s "- 1 } 

{%,Ki),«i,io):l<J<«il<iO')< r i}\{ u i.i> t 'i,i} 

■E{[i? x ] 

" " * ' 

[R* \mim- 2 ,k 1 k'2 

[-Rx]l)„ >1 D„ >2 ,M„ >1 M„ >2 [^x]tl„ j2 D„,3,tl„,2Mn 1 3 ' ' ' [^] V n _ rn V± _ j ,U„ Ui, i } ■ (H7) 

' " " ' 

[fii"]m„mi,fc„fci 

Denote the summation comprised of the second to fourth lines of (117) as ^4(«i,i, vi,i; {si}^ 1 ; {^}" = i)- 
We are going to show the limiting value of A(u^i, v\y, {si}^ 1 ; {^i}"=i) is 0(1)- Let 

s(Ki,« 3 ,i}"=2; Mr=i) = E {( rf «2,i ( u 2,i)) Sl (d« 3 ,i (^i)) 52 • • • K.Ki)) 8 - 1 } 

and 

c, ({ u j,i(j)^j,Ki)}i<i<n,i<Ki)<'' 3 ; *L ri ^=i) 
= E{[R X ] 

fl,lfl,2,«l,l«l,2 [-^x]t)l,2t)l,3,ttl,2ttl,3 ' ' ' [Rx]vi,r 1 V2,l,U ll r 1 U 2 ,l 

[^l]«n,lHn,a,«tl,lIltl,2[^l]l'n,2«i. 1 3,«lii,2ll.,3 ' ' ' ["^*] f n ,r„ «1, 1 ,"r.,7-„ «1 , 1 I ' 

That is, B({uj : i,Vj t i}; {si}) and C({uj t i^, Vj^}; {r^}) stand for the first and second expectations in 
^4.(^1,1,^1,1; {si}', {ri}), respectively. We have 

|*4(wi,i, v 1A ; {si}; {ri})\ 

< E \B({u jt i,v jtl y,{ 8i })\ ■ \C({u m) ,v m) };{r,})\ .(118) 

{u 3 ,i U ),v jtl(j y.l<j<n,l<l(j)<r j }\{u 1A ,v 1<1 } 

Due to Holder's inequality and bounded moments of dfc(m)'s, we have 

\B({u hU v hl }-{s t })\ < E{|K 2 , 1 h,))r i(n - 1) } 1/(f! - 1) ---E{|(4„, 1 (vi))| S "- l(n - 1) } 1/(n - 1) 

= 0(1)- (119) 
When x = cs, suppose that \B({uj : i,Vj : i}; {si})\ is less than a constant c\. Thus, by (118), 

<ci Yl \ c ({ u j,iu)' v j,m}'{ ri V\ • 

,i'i,iO):l<j<n,l<iO')< r i}\{ tt i.i, t 'i.i} 

Following the proof of Lemma 1 given in Appendix II, we have 

{uj,n ] ),v ] , HJ) :l<j<n,l<l(j)<r ] }\{u ltl ,v 1A } 
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which is equal to fi(RH A hr ") in the limit of K,M,N — > oo and if/iV — > /3 for any v± : i and ui 5 i. 
Thus, we have shown the limiting value of A(u\ t \,vij_\ {ri}) is O(l) when x = cs. 
When x = ca and condition 2) in Theorem 6 is true, we can show 

£ \C({u jm ,v jm };{n})\ = 0(1). (120) 

{ u 3,i(j), v j, i(j) :1<.?<»»,1<J0')< t V}\{ u i. i>"i,i} 

Thus, we can use (119) and (120) to demonstrate that limjv,A"=/9JV-»oo^4.(«i,i) v^i; {sj}; {r,}) = O(l). 
Consider the case that x = ca and condition 1) of Theorem 6 holds. Since (4(m)'s are non-negative 
random variables, B({uj^,Vj : i}; {s«}) is non-negative, and we suppose it is upper-bounded by constant 
C2- Then 

{uj,i U ),v ]MJ) :l<j<n,l<l(j)<r ] }\{u 1A ,v 1A } 
= C2^{[R c l + + n ]t)i,it)i,i,Ui,iUi,i}) 

which is equal to C2^{R r c i^ hr ") asymptotically. Thus, A(ui t i, vi t i; {sj}; {n}) is asymptotically equal 
to O(l) as well. 

It follows that, for Group 1 of noncrossing partitions of the ordered set {uj^ : 1 < j < n}, 
f i( Pl (R x )q 1 (D)--- Pn (R x )q n (D)) 
= lim (2M + l)- 1 ^ 1 ^6 nA Et^.Ki))'"} 

M,i\,K — ►oo 14 — * 14 — rf 

x ^ ai jri &i jSl ■ ■ ■a„ )rn 6 ri _i jSn _ 1 ^4(ui j i,ui j i;{s i };{ri}) 

ri,- ,r„ 

Si,— ,S n _i 

= 0. (121) 

Note that the second line of (121) is equal to b n tSn [i(D Sn ) = 0, while the third line is 0(1). Thus, 
in this case, p(pi(R x )qi(D) ■ ■ ■ p n (R x )q n (D)) = 0. 

Now, we consider Group 2 of noncrossing partitions of the ordered set {uj : i : 1 < j < n}. By 
Lemma 7, we can suppose that n t i and ut+i t i are partitioned in the same class. Under this circumstance, 
we can use Lemma 8 to show /j,(pi(R x )qi(D) ■ ■ ■ p n (R x )q n (D)) given by (112) and (113) can be written 
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as 

f i(p 1 (R x )q 1 {D)---p n (R x )q n (D)) 
= y^atrMK') lim (2M+1)- 1 K- 1 V ai, ri 6i, Sl • • • a„, r „6„, s „ x (122) 

^ — ' K,N,M^oo *■ — * s / 

x/iv-,/? ri ''"' r s t ;::. rt ;„ 1 ''"' r " ^^Z t 

£ £ E {fe(^)) Sl • • • (dk,^ (m t _i)) St - 2 (cZ fet (m t )) s *- 1+St (4 t+2 (m t+2 )) St+1 • • • 

K\{k t+1 } M\{m t+1 } 

(4j ( m l)) "}E{[i2 x 1 ] mim , 2i fc 1 fc 2 • • • [R^ 1 ]m t _ 1 m t ,k t _ 1 k t [R x t+1 ]rn t m t+ 2,k t k t+ 2 ' ' ' [R x "]m n m 1 ,k n k 1 } ■ 

Using similar arguments as in the discussion of Group 1, we can show the limiting sum behind J2 Tt a t,r t ^{R r K t ) 
in (122) is O(l) for x € {cs,ca}. Since at, rt iJ,(R r x ) = 0, for Group 2 of noncrossing partitions of 
the ordered set {uj : \ : 1 < j < n}, we have /j,(pi(R x )qi(D) ■ ■ ■ p n (R K )q n (D)) = 0. 

As we have shown both Groups 1 and 2 have contributions to /j,(pi(R x )qi(D) ■ ■ ■ p n (R x )q n (D)) equal 
to zero, the proof is completed. 
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